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Abstract 



The relation between the effective chiral Lagrangian whose parameters scale ac- 
cording to Brown and Rho scaling("BR scaling") and Landau Fermi-liquid theory 
for hadronic matter is discussed in order to make a basis to describe the fluctua- 
tions under the extreme condition relevant to neutron stars. It is suggested that 
BR scaling gives the background around which the fluctuations are weak. A simple 
model with BR-scaled parameters is constructed and reproduces the properties of 
the nuclear ground state at normal nuclear matter density successfully. It shows 
that the tree level in the model Lagrangian is enough to describe the fluctuations 
around BR-scaled background. The model Lagrangian is consistent thermodynam- 
ically and reproduces relativistic Landau Fermi-liquid properties. Such points are 
important for dealing with hadronic matter under extreme condition. On the other 
hand it is shown that the vector current obtained from the chiral Lagrangian is 
the same as that obtained from Landau-Migdal approach. We can determine the 
Landau parameter in terms of BR-scaled parameter. However these two approaches 
provide different results, when applied to the axial charge. The numerical difference 
is small. It shows that the axial response is not included properly in the Landau- 
Migdal approach. 
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1 Introduction 



Although QCD which deals with quarks and gluons is believed to be the 
fundamental theory for strong interactions, it is generally accepted that the 
appropriate theory at very low energy is the effective quantum field theory 
which incorporates the observed degrees of freedom in low-energy nuclear pro- 
cess, i.e., pions, nucleons and other low-mass hadrons [1-7]. The effective La- 
grangian in matter-free or dilute space is governed by QCD symmetries with 
its parameters to be determined from experiments in free space. The energy 
scale of the experiments in which one is interested determines which hadrons 
play an important role in the theory. For example, it is shown [7] that one can 
integrate out even the pions for the two nucleon systems at very low energy. 
According to the results of [7], the deuteron and low-energy nucleon-nucleon 
scattering properties can be described very accurately by an effective theory 
given in terms of the nuclonic degrees of freedom only with a cutoff around 
the natural scale of the theory which for low energy is the pion mass. 

Since there is currently growing interplay between the physics of hadrons and 
the physics of compact objects in astrophysics through the properties of hot 
and dense environments, we want to extend such successful strategy of effective 
field theories to a dense medium. First of all, our main goal is to understand 
the properties of dense hadronic medium which can be tested in various heavy 
ion collisions. Recent dilepton experiments (CERES and HEL10S-3) gave us 
very important information on the properties of hadrons in dense medium, 
which we will detail later. Understanding the properties of dense hadronic 
medium through heavy ion collision experiments is essential for understand- 
ing the properties of neutron stars which are believed to be formed in the 
center of supernovae at the time of explosions. Especially the determination 
of the maximum neutron star mass is one of the most important issues in astro- 
physics in explaining controversies between the observations and theoretical 
estimations, and it will give some hints for the detectibility of the gravita- 
tional wave detectors. A dense matter makes the probed energy scale larger 
than that in free space. Therefore we have to introduce more massive degrees 
of freedom which are usually vector meson and/or higher order operators in 
the nucleon fields. We must also consider a new energy scale, Fermi energy of 
nucleons in bound system. 

The standard strategy to attack the dense hadronic matter is to obtain the 
ground state of matter and compute the excitations around it, based on an 
effective Lagrangian whose parameters are obtained in free space. Although 
such an approach can give satisfactory results with a sufficient number of 
parameters, it is not obvious whether we can extend the results for the dif- 
ferent density. When the extension does not work, we usually face with very 
complicated loop diagrams which may lead to the impasse. 



7 



Another strategy is to start from the in-medium Lagrangian which is built 
on the reasonable assumptions, instead of deriving the hadronic matter prop- 
erties from the matter-free effective Lagrangian defined in free space. In this 
approach one regards its mean field solution as a solution of the Wilsonian 
effective action in which the high energy modes are integrated out into the 
coefficient. We can compare it with the well-known Landau Fermi-liquid the- 
ory, which works at low energy excitation in strongly correlated Fermi system. 
Landau Fermi-liquid theory is described by quasiparticles, which are the low 
energy excitations in Fermi liquid, and their interactions under the assumption 
of one-to-one correspondence between the quasiparticle in the liquid and the 
particle in a non-interacting gas. It is a fixed point theory [8,9] under Wilsonian 
renormalization to the Fermi surface [10] with A/kp — » if Bardeen-Cooper- 
Schrieffer(BCS) instability does not exist, where A is the cutoff of the theory 
relative to the Fermi surface and kp is Fermi momentum. Since the result after 
the repeated renormalization does not depend on hp, the argument for Fermi 
liquid holds as long as there is no phase transition. A famous example of such 
a Lagrangian is Walecka model. Its extension and justification were studied 
recently by Furnstahl et al. [11-13]. Their Lagrangian is constrainted by QCD 
symmetry and by Georgi's naturalness condition and the coupling constants 
in it are tuned in order to describe nuclear ground state. Bulk properties of 
nuclei are described by it very successfully. But it is somewhat unclear to know 
how to approach the fluctuation on the ground state. 

In this review we will approach the nuclear matter in a different way. We wish 
to apply the strategy of the effective field theory to the in-medium theory. 
We assume that the in-medium effective Lagrangian has the same structure 
as in free space according to the symmetry constraint of the fundamental the- 
ory QCD but that its parameters are modified in medium. It means that the 
effect of embedding a hadron in matter appears mainly in the change of the 
"vacuum," i.e., quark and gluon condensate in QCD variables and the param- 
eters in an effective theory. In this strategy the density- dependent parameters 
include many-body correlations. 

Brown-Rho(BR) scaling [14] is one specific way to define such in-medium 
parameters. Brown-Rho scaling is the scaling of the dynamically generated 
masses of hadrons which consist of chiral quarks, i.e., u and d quarks. Brown 
and Rho phrased the scaling with the large N c Lagrangian, i.e., Skyrmion, 
under the assumption that the chiral symmetry and the scale symmetry of 
QCD are relevant. Their Lagrangian is implemented with scale anomaly of 
QCD, too. The masses and pion decay constant of this QCD effective theory 
scale universally: 



(1) 
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The star represents in-medium quantities here, v is vector meson degree of 
freedom and s an isoscalar scalar meson which has a mass ~ 500 MeV in 
nuclear matter. M represents a free nucleon mass and M* a scaled nucleon 
mass which is somewhat different from the Landau effective mass discussed in 
Section 6.1. 

It is known that BR scaling describes that the light-quark vector meson prop- 
erty in the extreme condition very successfully. One can make the extreme 
condition through relativistic heavy ion collisions. Specially the dileptons pro- 
vide a good probe of the earlier dense and hot stage of relativistic heavy 
ion collision since the interaction of leptons is not subject to the strong in- 
teractions of the final state. CERES (Cherenkov Ring Electron Spectrometer) 
collaboration observed in heavy ion collision (S+Au) that the dilepton pro- 
duction with invariant masses from 250 MeV to about 500 MeV is enhanced 
much more than the predicted from the superposition of pp collision [15]. And 
HELIOS-3(CERN Super Proton Synchrotron detector) [16] also observed the 
dilepton enhancement in S+W. It is shown by Li, Ko, and Brown [17] that a 
chiral Lagrangian with BR-scaled meson masses describes most economically 
and beautifully the enhancement, which is found to come from the dropping 
vector meson masses in dense matter. Furthermore the excitation into the 
kaonic direction above the given ground state seems to describe the properties 
of kaons in medium [18,19] with scaled parameters [20]. 

Since BR scaling gives the universal scaling mass relation among hadrons, it 
must work for the baryon properties in medium. In recent works [21-25] it 
has been discussed how the BR scaling, which is applied to meson properties 
successfully, can be applied to the baryon property in dense medium and 
how it can be extrapolated to a hadronic matter under extreme conditions 
from the known normal nuclear matter. The construction of such bridges will 
be necessary to understand the various phenomena in relativistic heavy ion 
collisions and in compact stars in the universe. For those purposes we develop 
arguments for mapping the effective chiral Lagrangian whose parameters are 
governed by BR scaling to Landau Fermi-liquid theory. We will relate the 
meson mass scaling with baryon mass scaling and also relate matter properties 
with BR scaling parameter, which implies the vacuum structure characterized 
by quark condensate, via Landau parameter. Fermi liquid theory may work up 
to chiral phase transition, though the relevant degrees of freedom are changed 
to quasiquarks from quasihadrons. 

In this review we approach the aim in two ways. The first is to obtain the 
ground state with BR scaling. How the Fermi surface is obtained in BR scaling 
framework is not yet understood clearly. So people usually assume that the 
ground state of hadronic matter is determined by the conventional matter 
from a standard many-body theory and use density- dependent effective chiral 
Lagrangians to compute mesonic fluctuations above the ground state. Though 
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various fluctuation phenomena can be described successfully in this way, such a 
treatment gives no constraints for consistency between the excitations and the 
ground state. Though BR scaling is applied very successfully to describe meson 
properties in medium [17,26], the way to deal with the matter properties is 
disconnected from BR scaling. Such a procedure is not satisfactory in going to 
the higher density region from the normal nuclear matter density. We can take 
the kaon condensation in neutron star as an example. In dealing with it, we 
come across a change of the ground state from that of a non-strange matter to 
a strange matter. The works up to date [27,28] treated KN interaction and the 
ground state separately. This is not satisfactory. Ground state properties might 
effect the condensation. For example, the effect of the four-Fermi interactions 
which play an important role in determining ground state, suppresses a pion 
condensation [29]. So we must deal with the whole bulk involving the ground 
state and excitations on top of it on the same footing. We assume that the 
ground state is given by the same effective Lagrangian that is supposed to 
include higher order corrections as the mean field of the BR-scaled chiral 
Lagrangian. We want to make an initial step for dealing with the ground 
state and the fluctuation on top of it on the same basis. We bridge these 
two properties by constructing a simple model whose parameters scale in the 
manner of BR and which describes nuclear matter properties well. It is shown 
that the model can be mapped to Landau Fermi-liquid theory. 

The next step to achieve our aim is to identify the parameters of the BR-scaled 
effective Lagrangian with the fixed point quantities in Landau Fermi-liquid 
theory, given a hadronic matter with a Fermi surface. With this identification, 
certain mean field quantities of heavy meson (e.g. p, a;) can be related to BR 
scaling through the Landau parameters. We show how such arguments work 
for the electromagnetic currents in nuclear matter. We can link a set of BR- 
scaled parameters at nuclear matter density with the orbital gyromagnetic 
ratio in terms of the Landau parameter which comes from integrated-out 
isoscalar vector degree of freedom in the effective Lagrangian. Then we will 
try to derive the corresponding formulas for the axial current in a similar way. 

This review is organized as follows. In Section 2 Landau Fermi-liquid theory 
and its interpretation in terms of renormalization group language are summa- 
rized briefly. And we show how thermodynamic observables are related to the 
relativistic Landau parameters. In Section 3 the strategy of an effective field 
theory and how it is applied to a dense matter are explained. We proceed to 
discuss how nuclear matter described by an in-medium effective Lagrangian 
can be identified with Landau Fermi liquid. The model of Furnstahl, Serot 
and Tang [11] (referred to FTS1) which imposes the symmetries of QCD is 
examined as an example of the application of a general strategy of an effective 
chiral theory to a medium. In Section 4 Brown-Rho scaling is derived with 
QCD-oriented effective Lagrangian. A model where BR scaling governs the 
parameters of a chiral Lagrangian and determines the background at a given 
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density is constructed in Section 5 in order to describe in weak coupling the 
same physics as FTS1 which has strong coupling in the form of a large anoma- 
lous dimension of a dilaton. It describes well normal nuclear matter properties 
and has thermodynamic consistency and Fermi-liquid structure needed for 
the extrapolation to higher density region. In Section 6 vector current and 
axial charge transition matrix elements for a nucleon above the given Fermi 
sea in Landau-Migdal theory and in chiral effective Lagrangian are calculated 
and compared. A summary and comments on some unsolved problems are 
given in Section 7 and Section 8 respectively. Appendix A shows how sen- 
sitive the equation of state is to the many-body correlation parameters for 
p > pq. Appendix B shows how to compute relativistically the pionic contri- 
bution to Landau parameter F\ by means of Fierz transformation. And the 
vector-mesonic contribution to F\ and to electromagnetic convection current 
is calculated in Appendix C relativistically with random phase approximation. 



2 Landau Fermi-liquid theory 

We discuss briefly Landau Fermi-liquid theory in this section, before presenting 
the relation between the chiral effective theory for nuclear matter and Landau 
Fermi-liquid theory. A mini-primer on Landau Fermi-liquid theory is given in 
the first subsection to define the quantities involved. In the second subsection, 
it is discussed that Landau Fermi-liquid theory is considered as an effective 
theory and is shown to be a fixed point theory in renormalization group(RG) 
language. And we will show how the thermodynamic quantities are related to 
relativistic Landau parameters in Section 2.3. 



2.1 A mini-primer 

Landau's Fermi-liquid theory is a semi-phenomenological approach to strongly 
interacting normal Fermi systems at small excitation energies. The elementary 
excitations of the Fermi-liquid, which correspond to single particle degrees of 
freedom of the Fermi g called quasiparticles in Landau Fermi-liquid 

theory. It is assumed that a one-to-one correspondence exists between the 
low-energy excitations of the Fermi liquid near Fermi surface, i.e., quasipar- 
ticles, and those of a non-interacting Fermi gas. A quasiparticle state of the 
interacting liquid is obtained by turning on the interaction adiabatically at the 
corresponding state of non-interacting Fermi gas. The quasiparticle properties, 
e.g. the mass, in general differ from those of free particles due to interaction 
effects. In addition there is a residual quasiparticle interaction, which is pa- 
rameterized in terms of the so called Landau parameters. 
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The adiabatic process described above is possible in the vicinity of the Fermi 
surface only. Let us see the Fermi liquid at T = 0. Since Pauli exclusion princi- 
ple makes the states below the Fermi surface filled, quasiparticle with energy 
e loses energy less than e — e F when colliding the background particles. It 
means that the quasiparticles which can interact with the quasiparticle are 
those with an energy within \e — Sf\ of the Fermi surface. And the final state 
momenta are also restricted by e' < e. Pauli exclusion principle and the cor- 
responding rarity of final states make the quasiparticle life time proportional 
to \e — e F \~ 2 at T = case. 

Fermi-liquid theory is a prototype effective theory, which works because there 
is a separation of scales. The theory is applicable to low-energy phenomena, 
while the parameters of the theory are determined by interactions at higher 
energies. The separation of scales is due to the Pauli principle and the finite 
range of the interaction. Pauli principle makes the low energy quasiparticle 
physics possible near the Fermi surface and the finite range of interaction 
makes a few quasiparticles around Fermi surface, who appear by the small 
change of the energy in low energy physics, form a gas. Fermi-liquid theory 
has proven very useful [30] for describing the properties of e.g. liquid 3 He and 
provides a theoretical foundation for the nuclear shell model [31] as well as 
nuclear dynamics of low-energy excitations [32,33]. 

The interaction between two quasiparticles p 1 and p 2 at the Fermi surface of 
symmetric nuclear matter can be written in terms of a few spin and isospin 
invariants [34] 



fp 1 a 1 r 1 ,p 2 



N(0) 



F(cos 9 12 ) + F'(cos 9 12 )r 1 • r 2 + G(cos 9 12 )cr 1 ■ er 2 

kp 



q 2 

+ G'(cos9 12 )cr 1 • <t 2 t x ■ r 2 + -^H (cos 9 12 )S 12 (q 



+ |!#'( cos 12 )S 12 (g) Tl -T 2 

Kp 



(2) 



where 9yi is the angle between p x and p 2 and N(0) = * s ^ ne density 

of states at the Fermi surface. In this review natural units where h — 1 are 
used. The spin and isospin degeneracy factor 7 is equal to 4 in symmetric 
nuclear matter. Furthermore, q = p x — p 2 and 

S 12 (q) = 3<Ti • qcr 2 ■ q - a x • cr 2 , (3) 



where q = q/\q\- The tensor interactions H and H' turn out to be important 
for the axial charge [34] The functions F, F', . . . are expanded in Legendre 
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polynomials, 

F(cos0 12 ) = ^F,P,(cos0 12 ), (4) 

t 

with analogous expansion for the spin- and isospin-dependent interactions. 
The energy of a quasiparticle with momentum p = \p\, spin a and isospin r is 
denoted by e P)(TjT and the corresponding quasiparticle number distribution by 
n Pt a,r- From now on the spin and isospin indices a and r will be omitted from 
the formulas to avoid overcrowding, except where needed to avoid ambiguities. 



n p (r,t)=n° p (e° p ) + 5np(r,t) (5) 
j^fpp>6n p ,(r,t) (6) 

where 8n p (r,t) is the long wave length excitations from the ground state 
n,p(e°) in the vicinity of Fermi surface. The space and time dependence of the 
quantities will also be omitted, e.g., e = e(r,t). The Landau effective mass 
and velocity of a quasiparticle on the Fermi surface is defined by 



P = Kp L, 



dp 

We must note that the total current is not 



/ ''" ~v* F . (7) 



JlocQP = E / 77^ V F 6n P ( 8 ) 



d 3 p 

in Landau Fermi-liquid theory. The quasiparticle distribution n p (r,t) obeys 



The key assumption of Landau Fermi-liquid kinetic theory is that e p (r, t) plays 
the role of the quasiparticle Hamiltonian; 



f t =-VrS p (10) 

dv 

The kinetic equation is 
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+ V r n p ■ V p e p - V p n p ■ V r e p = I[n p ] (12) 

where I[n p '] is an internal collision integral which represents the sudden change 
of quasiparticle momenta. We consider the system without external forces. 
Integrating over p, that cancels the effect of the internal collision under the 
assumption that the quasiparticle energy, momentum and number are locally 
conserved, (12) becomes (to order 5n p (r,t)), by integration by part, 

/d v OSti OS o f d V 

= o (1 

where 5p is the total particle (or quasiparticle) density fluctuation. We can 
easily see that the conserved current is 



<T,T 



1 dp 



E/tS^H" (14) 



d 3 p 

(2k)* 
d 3 p 
(2^ 

with the excitation from the local equilibrium n°(e(r, t)); 

np(r,t)=n p (e(r,t)) + 5n l p oc (r,t). (15) 
Comparing (15) with (5) and (6), we obtain 



The modification Ji oc qp — J can be interpreted as the effect of the return 
flow of the surrounding matter to the localized wave packet carrying J and is 
called back-flow current. 



2.2 Renormlization group approach 



Landau Fermi-liquid theory was based on Landau's reasonable intuition. Af- 
ter his work, his theory was derived and justified microscopically [35]. Recent 
development of Wilsonian RG method in medium [10] provides us a new un- 
derstanding of Fermi liquid theory. Landau Fermi-liquid theory is a fixed point 
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Fig. 1. Our interest lies in the spherical shell which is of the thickness of A on either 
side of Fermi sphere. In RG transformation we eliminate the mode within the two 
thick-lined shells whose thickness is dA. 

theory described by marginal coupling, i.e. m* L [8]. In this section we will 
review the RG arguments for Fermi liquid theory. 

To do this we consider a nonrelativistic system of spinless fermions whose 
Fermi surface is spherical characterized by kp for simplicity. Then non-interacting 
one particle Hamiltonian near Fermi surface is 

H = ^— — — —k F = v F k (17) 
2m 2m m 

with k = \K\ — kp. The free fermion field action 

S = J rjj(ukQ)(iu - vk)tp(u}kQ) (18) 

A 

in momentum space where 



r = r _dn_ f dk_ J du_ 

I ■ J (2nyJ A (2n)_l (2k) ' 1 J 

■0 and ip are Grassmannian eigenvalue with fermion operator ^; ^\ip) = 
and (^l^t = Note that a shell of thickness A on either side of the 

Fermi surface is taken for low energy physics as seen in Fig. 1. Pauli exclusion 
principle lets only small deviations from the Fermi surface, not from the origin, 
be important in low energy physics of fermion matter. This defines the starting 
point of an in-medium renormalization group procedure. 

The first step for the renormalization is decimation; to integrate out the high 
energy mode whose momentum is larger than A/s and to reduce the cutoff 
from A to A/s. For example, the free action (18) becomes 
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A/s oo 
—A/s — oo 

Next we rescale the momenta in order to compare the old and the new; 



(u,k) -> (suj,sk). (21) 
The last step is to absorb the uninteresting multiplicative constant by 

if) -> S - 3/2 i{j. (22) 

The RG transformation consists of these three steps. After such RG trans- 
formation, the free action (18) returns to the old. When a coupling is turned 
on, the coupling is called relevant if it increases after RG transformation. If it 
decreases, it is called irrelevant and if it remains fixed like the free action, it 
is called marginal. 

Now let us turn on four-Fermi interaction. The appropriate action is 



S = J ip[iw — Vpk]t/j + 5 pi* J ipip 

A A 

+ 2!2~! / M(4 ' 3 ' 2 ' 1 )^(4)^(3)^(2)^(1) (23) 
'a 4 

where (i) represents (u^, fcj, Qj) and 



/-n/p^/p^/^IA-W). (24) 

A 4 1 —A — oo 

with a cutoff function for k±, 0(A — |&4|), which is needed to make all the 
momenta lie within the band of width 2A around the Fermi surface. Here 
v F = kp/m* and 



m * = — _ (25) 

is the effective mass of the nucleon which will be equal to the Landau mass m* L 
as will be elaborated on later. S(cj, k) is self-energy and l/Z = 1 The 
effective mass arises because the eliminated mode contributes to iujipip and 
kipip differently. By defining tp' = s~ 3//2 Z -1 / 2 -^ we fix the coefficient of iwijjij) 
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Fig. 2. Momentum conservation and the physics near Fermi surface force the angle 
between K% and K4 to be the same as that between K\ and K2 unless K2 = —K\- 



and define the effective mass. The term with 5/j,* is a counter term added 
to assure that the Fermi momentum is fixed (that is, the density is fixed). 
What this term does is to cancel loop contributions involving the four-Fermi 
interaction to the nucleon self-energy (i.e., the tadpole) which contributes 
marginally so that the v* F is at the fixed point. Since other contributions 
which can be written as —Htpip are irrelevant, m moves to m* in earlier stage 
of renormalization but becomes a fixed point characterized by some m*. This 
means that the counter term essentially assures that the effective mass m* be 
at the fixed point. Without this procedure, the term quadratic in the fermion 
field would be "relevant" and hence would be unnatural [8]. 

Let us see the quartic coupling u at tree level. The cutoff function 0(A — |fe 4 |) 
in (24) makes the coupling depend on angles on the Fermi surface. Since all 
the momenta are on the thin spherical shell near fop, momentum conservation 
makes the new quartic coupling u'(uj', k', Q') decay after the RG transforma- 
tion at tree level except for the two cases. One is that K% and K4 are the 
rotation of Ki and K 2 around the K ± + K 2 as seen in Fig. 2. In that case, the 
opening angle cos -1 ^! • U2) is fixed where ui is a unit vector in the direction 
of Ki. The other, so-called BCS coupling, is that u\ = —U2 and = — U4. 
These cases are represented by functions; 



where 6ij = Ui ■ iij and $ is the angle between the planes containing (Ki, K 2 ) 
and (-K3, K4) respectively. T and V are marginal at tree level. 

Then the next question is how T and V evolve at one loop level. Figure 3 shows 
the one loop diagrams for the evolution. Integrating out the momentum shell 
of thickness dA at k = ±A and sending A//cp — > 0, all the diagrams in Fig. 3 
do not contribute to T . So T is left marginal to the one loop order. In the case 
of V, the third diagram in Fig. 3 makes a flow. When we expand V in terms 



«(012 
«(013 



034)= ^(012, #) 

M = v(M 



(26) 
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Fig. 3. One loop diagrams for the renormalization of the marginal quartic couplings 
T and V. 

of angular momentum eigenfunction, V becomes irrelevant only if all the V;'s 
are repulsive. However, if any Vi is attractive, it becomes relevant and causes 
BCS instability. We call it BCS channel. Since Landau theory assures that 
there is no phase transition for one-to-one correspondence between particles 
and quasiparticles, BCS channel destabilizes in Landau Fermi-liquid theory. 

If we divide the angular part of the shell integration in the action (18) into the 
cells of size Af2 ~ A/k F , the shell is split into ~ k F /A cells. By analogy with 
large N theory each cell corresponds to one species, i.e., N ~ k F /A. l/N ex- 
pansion tells that all higher loop corrections except for the bubbles of the first 
diagram in Fig. 3 are down by powers of l/N = A/k F . And the contributions 
of bubbles of the first diagram survive but are irrelevant to the (3 function. 
Only one loop contributions to f3 function survive in the large N limit. So 
the four-Fermi interactions in the phonon channel T are also at the fixed 
points in addition to the Fermi surface fixed point with the effective mass m*. 
Note that only forward scattering T{& = 0) is important for responses to soft 
probes, since nonforward amplitudes in loop calculations are also suppressed 
by l/N. Six- Fermi and higher- Fermi interactions are irrelevant and contribute 
at most screening of the fixed-point constants. The Landau parameter F can 
be identified with forward scattering 7 ' — 0) easily [8]. We arrive at 

the Fermi-liquid fixed point theory in the absence of BCS interactions. 

Chen, Frohlich, and Seifert [9] obtain the same result in the l/N expansion 
where their N is taken to be N ~ A with 1/A being the width of the effective 
wave vector space around the Fermi sea which can be considered as the ratio of 
the microscopic scale to the mesoscopic scale. More specifically if one rescales 
the four-Fermi interaction such that one defines the dimensionless constant 
g ~ u /k F where u is the leading term (i.e., constant term) in the Taylor 
series of the quantity u in (23), then the fermion wave function renormalization 
Z, the Fermi velocity v F and the constant g are found not to flow up to order 
0(g 2 /N). Thus in the large N limit, the system flows to Landau fixed point 
theory to all orders of loop corrections. This result is correct provided there 
are no long-range interactions and if the BCS channel is turned off. 
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2.3 Relations for relativistic Fermi-liquid 



In this section we briefly summarize the relations of physical properties of 
the relativistic Landau Fermi-liquid. The extension of Landau Fermi-liquid 
theory to relativistic region for high density matter is found in [36]. It should 
be pointed out that one can use all the standard Landau Fermi-liquid relations 
established below in our calculations once fixed point quantities are identified 
in the chiral Lagrangian. 



2.3.1 Compression modulus and F Q 
The density of states at Fermi surface is 




Note that m* L = y m 2 N + kp for relativistic non- interacting Fermi gas. m* L 
defined by (7) includes the kinetic energy for relativistic Fermi liquid. The 
chemical potential is defined by 



d£ , , 

^ Tp =s F (28) 

where S represents the energy per volume. Using (27) and (28) one can derive 
the relativistic relation between compression modulus which represents the 
change of volume with pressure and F in the same way as the nonrelativistic 
one: 



djjj ( f \ 
K = 9p d^ = 9p lh V° p + J dT ' f pp' n p'J 



dp dp 



3k 



m. 



1 + F ) (29) 



Here p is the baryon number density and n p is the Fermi distribution function 
for the state of momentum p. 



2.3.2 Landau effective mass 

In deriving the Landau mass formula, we shall compare the rest frame and the 
boosted frame with very small velocity u and check Lorentz symmetry. Since 



19 



u = I u I is small we neglect the order of u 2 in the process of the derivation. 
(Note that 7„ = (1 — w 2 ) -1 / 2 1 in that order. ) 

Let us first derive the relativistic relation between the Landau effective mass 
and the velocity dependence of the quasiparticle interaction. When we add a 
particle (or equivalently quasiparticle) of momentum p in the rest frame to 
the system, the energy and the momentum of the system increases by p and 
s p (0) respectively. In a moving frame with velocity — u, the momentum and 
energy increase by 



p' = p - u(u ■ p)(l - 7u ) + £p (0)u7, 
^p + e p (0)u 
e p ,(u) = (e p (0) + p-u)^ u 
^E p (0) + pu. 



u 



(30) 



(31) 



From e p (u) 



(30) and (31), we have 



e p (u) 



£p- £p (o) u (0) + p ■ u = e p (0) - e p (0)u ■ V p £p(0) + p ■ u. 



(32) 



and 




(33) 



Since n p i{u) = n p (0), we can obtain 



n p (u) w n P - £pU (0) = n p (0) - e p (0)u ■ V p n p (0) 



(34) 



using (30) and (31). Then (33) becomes 




(35) 



Comparing it with (32) 




(36) 



In the ground state 



Vp7lp(0) 



S(e p 




(37) 
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gives 



P = £p ~dp + J ' dr6 ( £ P ~ v) £ pfppP ' P-fip- ( 38 ) 
The chemical potential \x is e PF . Thus (36) becomes on the Fermi surface 



(de p \ IPph { on\ 

»{-^) PF =PF -^T (39) 

Using (27), one obtains 



= ! + (40) 

This is the relativistically extended relation of the famous Landau mass for- 
mula 

^ = 1 + ^, (41) 
M 3 K J 



that follows from Galilean invariance in the same way [30]. 



2.3.3 First sound velocity 

The first sound is a density oscillation mode under the circumstances where 
there are many quasiparticle collisions during the time of interest. The suffi- 
ciently often collisions produce the required local equilibrium in the time scale 
of the period of motion. 

When first sound wave gives small change in density of a static homogeneous 
relativistic fluid in a comoving frame without changing entropy, the continuity 
equation requires 



J^p = -pV • v. (42) 
Under the condition of the fixed entropy 



kds = —(pSv + Se) 



p 2 T 



{pSe - (p + S)5p} = 



(43) 
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with the entropy per particle ks, the volume per particle v = 1/p, and the 
energy per particle e = S/p, the relatvistic equation of motion becomes 



dv 
~dt 



p + S 



1 



VSp = 



P 5S 



5p 



V5p. 



(44) 



Applying (44) to (42), we obtain 



dt 2 




V 2 5p. 



(45) 



Thus the first sound velocity of the relativistic Landau Fermi-liquid is 



2= dp_ 

1_ dS 

dp d ( p2 d£/ P \ 
dS dp \ dp J 
K 
~9p 

k 2 F 1 + F 
3^1 + ^/3 



from the above results. 



3 Chiral effective Lagrangian for nuclei 

In this section we come to hadron/nuclear phenomena. We simply explain 
the basic theory of strong interaction, QCD, and how successfully low energy 
phenomena can be described by its effective theory. Then it is discussed how 
the effective theory can be applied to the nuclear /hadronic matter ground 
state. FTS1 model will be studied as an example. 

3.1 QCD: basis of strong interactions 

Quantum Chromodynamics is believed to be the fundamental theory for the 
strong interaction. It is a nonabelian gauge description of the strong inter- 
action and its building blocks are quarks and gluons. The QCD Lagrangian 
is 
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C QCD = -^Tr G llv G' iV + q(iql - gG- m q )q (47) 

with q = (u,d,...) T and m q =diag[m M , rrid, ■■■} in flavor space. In this section 
Tr represents the trace over flavor indices. The gauge field strength tensor is 



Gn V = d^Gy - d v Gfj, - ig[G„ G v \ (48) 

where G^ = G a ^. \ a is the well-known Gell-Mann's 3x3 traceless hermitian 
matrix with color indices a = l,2,...,8. 

Quantum Chromodynamics conserves its relevant symmetries; Lorentz sym- 
metry, parity, charge conjugation, time reversal and color SU(3) gauge symme- 
try. In addition, it has approximate symmetries. Isospin symmetry is a good 
one because of small difference between u and d masses. 

There is another well-known and important QCD approximate symmetry, chi- 
ral symmetry. It is based on the fact that the mass of light quarks u and d 
are much smaller than chiral symmetry breaking scale A x ~ 471"/^ ~ 1 GeV. 
We can consider the massless limit, i.e. m q — > when dealing with light quark 
physics. It means that chiral symmetry becomes relevant. With projection 
operator 



P± = (49) 
we can decompose the spinor into the eigenstates of helicity <x • k, 

akq = ±q±. (50) 
with q± = P±q. Without quark mass term, q + and g_ are decoupled each other 

C QCD - —TrG^G^ + q + {icf - gG)q + - - gG)q^ (51) 

So the separate transformations 

q± - e~^q ± (52) 

and 



q± - e~ ie± - t q± (53) 
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with real parameters e± and 6 p m leave the massless QCD Lagrangian invari- 
ant, where t is the general SU(2) generators. We can construct SU(2) and 
U(l) vector transformations as 



q -> e- i0y t q, q -> e~^g (54) 
and axial ones as 



q -> e'^^q, q -> e- ieA75 g (55) 
with = 0y ± #A and e± = ey ± e^- The corresponding currents are 



V^q^tq, V? = itfq, (56) 
A» = qYl 5 tq, A? = <rf»>y 5 q. (57) 



Although massless QCD Lagrangian is invariant under these transformations, 
its U(l) axial current of QCD is not conserved due to anomaly even if the 
quark mass m q — > 0; 



d^ s = ^Tr G^G*" + 2tqm ql5 q (58) 

with = e^G aP . So the massless QCD has chiral symmetry (SU(2) fl xSU(2) L ) 
and fermion number symmetry (U(l)y). Chiral symmetry does not appear in 
hadron spectrum. So it is generally assumed that the chiral symmetry is dy- 
namically broken into SU(2)y in hadron physics and the light pseudoscalar 
mesons are regarded as Goldstone bosons. 

Since the massless QCD Lagrangian contains no dimensional parameters, 
massless QCD action is invariant under scale transformation; 



x = ax, q(x') = a- 3/2 q(x), G^x') = a^G^x). (59) 

Thus massless QCD seems to have another approximate symmetry, scale sym- 
metry. However, the renormalization prescriptions which have the running cou- 
pling constant break the scale symmetry. It means that we need to introduce 
a dimensional scale in order to specify the value of running coupling constant. 
Such specification of a scale breaks the scale invariance of QCD; 



8»D» = 6£= ^Tr G, U G^ + q(l + lq )m q q (60) 
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with the anomalous dimension of quarks 7 g =diag[7„, 7^, ...] even if quark 
masses are zero. is dilatation current, 8% is the improved energy-momentum 
tensor [37] and 



with the scale of renormalization fi. Contrary to the axial anomaly which 
includes only one-loop contributions, scale anomaly includes multi-loop con- 
tributions in /3. Scale anomaly which plays a major role in breaking scale 
invariance of light quark QCD gives a basis of BR scaling argument which will 
appear in Section 4.2. 



3. 2 Effective field theories 

We know that the particles which appear in nuclear physics, e. g. pions, nu- 
cleons, etc., are not elementary particles. In the framework of QCD, they are 
known to consist of quarks and gluons which give non-perturbative contri- 
butions in low energy processes. In addition, QCD is a part of the Standard 
Model which is also an effective theory, not a fundamental theory. 

In order to deal with low energy hadrons, one can construct an effective field 
theory that is appropriate to the probed energy scale Q. The relevant degrees 
of freedom in such effective field theory are the low-energy particles which 
appear at the energy scale of the observed experiments. The particles with 
energies higher than the relevant energy scale are integrated out and absorbed 
into the couplings among the relevant degrees of freedom; 



where 77, fj, and j are the external sources of elementary fermions g(quarks), 
anti-fermions g(anti-quarks), and bosons g (gluons) respectively in fundamen- 
tal theory(QCD) and B, B, and M represent fermions (baryons), anti-fermions(anti- 
baryons), and bosons (mesons) respectively in effective field theory. Although 
one can in principle calculate the couplings in effective theory for strong in- 
teractions from the fundamental QCD from equation (62), it seems an im- 
possible mission because we do not know how to perform such a highly non- 
perturbative calculation. However, we do not have to deal with, nor to know 
exactly, the fundamental theory. 




(61) 




(62) 
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The strategy to build an effective theory is simple. It consists of writing down 
the most general Lagrangian which conserves all the relevant symmetries that 
figure at a given energy scale, and satisfy the basic principles(e.g. quantum 
mechanics, cluster decomposition) of the theory. Then any theory under the 
same constraints looks like his/hers at sufficiently low energy scale though 
he/she cannot insist that the right theory necessarily leads only to his/hers. 
(This is called "folk theorem" by Weinberg [38] though he used it to explain 
the usefulness of quantum field theory.) By imposing the relevant symmetries 
of the QCD or of a more fundamental theory, e.g., chiral symmetry, and the 
basic principles we can build an effective field theory for low-energy nuclear 
processes where the composite hadrons are taken as elementary. 

The number of terms in the effective Lagrangian which are consistent with 
fundamental or assumed symmetries may be infinite, but one can manage to 
describe the probed physical processes by expanding the terms in Q/A, where 
A is s suitable cut-off scale and Q the scale probed (Q « A). So only a few 
leading terms are relevant for low-energy processes. And the couplings can be 
determined from available experiments. 

Such effective field theories with chiral symmetry so designed work very well 
in matter-free or dilute space. They are used to calculate soft-pion processes 
firstly [1], and recently extended to the process in light nuclei though nucleons 
are not soft [2-7]. Park et a/.'s work for two-nucleon systems at very low 
energies [7] is one of the best and newest examples of the success. Setting the 
cutoff near one pion mass, they integrate out all mesonic degrees of freedom, 
even the pions. Since the deuteron bound state is dilute, the parameters of the 
theory can be determined by free space experiments. The results in [7] confirm 
that the strategy of the effective field theory works remarkably well. When 
the pion field is included in addition, it provides a new degree of freedom and 
improves the theory even further allowing one to go higher in energy scale [7] . 

But in heavier nuclei, the energy scales of the system will be higher since 
the interactions between nucleons in such systems sample all length scales 
and hence other degrees of freedom than nucleonic and pionic need be intro- 
duced. It means that the irrelevant and neglected terms in the computation 
for the light nuclei become more relevant. We need to consider more and more 
terms as the density of the system becomes higher and higher. The power of 
effective field theory whose parameters are determined from the matter-free 
experiments gets weakened in dense system. How can we proceed as density 
increases beyond the ordinary matter density for which there are practically 
no experimental data? 
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3.3 Effective Lagrangian in medium and Landau theory 



Recently Lynn made a progress to give a good hint to answer the above 
questions [39]. He proposed that the ground-state matter is "chiral liquid" 
which arises as a non-topological soliton. Fluctuation around this ground state 
should give an accurate description of the observables that we are dealing with. 
We shall here extend this argument further and make contact with Landau's 
Fermi-liquid theory of nuclear matter. This will allow us to understand the 
nuclear /hadron matter description in terms of chiral Lagrangians and Fermi- 
liquid fixed point theory thereby giving a unified picture of ordinary nuclear 
matter and extreme state of matter probed in heavy-ion collisions. This is the 
attempt to connect the physics of the two vastly different regimes. 

The basic assumption we start with is that the chiral liquid arises from a 
quantum effective action resulting from integrating out the degrees of freedom 
lying above the chiral scale A x ~ 47r/ 7r ~ 1 GeV. 



where the subscript <(>) represents the sector u < A x (u > A x ) of the given 
set of fields 0. As explained in Sec. 3.2, 



is the sum of all possible terms consistent with symmetries of QCD. This 
corresponds to the first stage of "decimation" [9] in our scheme. The mean 
field solution of this action is then supposed to yield the ground state of nuclear 
matter with the Fermi surface characterized by the Fermi momentum kp. The 
effective Lagrangian was given in terms of the baryon, pion, quarkonium scalar 
and vector fields. The gluonium scalars are integrated out. Instead of treating 
the scalar and vector fields explicitly, we will consider here integrating them 
out further from the effective Lagrangian. This will lead to four-Fermi, six- 
Fermi, etc. , interactions in the Lagrangian with various powers of derivatives 
acting on the Fermi field. The resulting effective Lagrangian will then consist 
of the baryons and pions coupled bilinearly in the baryon field and four-Fermi 
and higher- Fermi interactions with various powers of derivatives, all consistent 
with chiral symmetry. A minimum version of such Lagrangian in mean field 
can be shown to lead to the original (naive) Walecka model [5]. 

The next step is to decimate successively the degrees of freedom present in the 
excitations with the scale E < A x [9]. To do this, we consider excitations near 
the Fermi surface, which we take to be spherical for convenience characterized 




(63) 



s x = Z 9iOi 



(64) 
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by kp. First we integrate out the excitations with momentum p > ±A (where 
p = \p\ and A < A x ) measured relative to the Fermi surface (corresponding to 
the particle-hole excitations with momentum greater than 2A). We are thus 
restricting ourselves to the physics of excitations whose momenta lie below 
2A as in Section 2.2. Leaving out the pion for the moment and formulated 
non-relativistically, the appropriate action to consider can be written in a 
simplified and schematic form as Eq. (23). R 

In nuclear matter, the spin and isospin degrees of freedom need to be taken into 
account into the four-Fermi interaction. All these can be written symbolically 
in the action (23). The function u in the four- Fermi interaction term therefore 
contains spin and isospin factors as well as space dependence that takes into 
account non-locality and derivatives. For simplicity we will consider it to be a 
constant depending in general on spin and isospin factors. Non-constant terms 
will be "irrelevant." In our discussions, the BCS channel that corresponds to 
a particle-particle channel does not figure and hence will not be considered 
explicitly. 

The successive mode elimination 

J [d^e* 5 *^ = J e ^<>^<), (65) 

which satisfies RG equation explained in Sec. 2.2, will give 



S* = £#0?- (66) 

i 

The starred coupling constant g* and operators O* depend on density through 
s and have the same structure of (64). 1(h) represents the components p < 
A/s(p > A/s). The upshot of the analysis in Section 2.2 is that the resulting 
theory is the Fermi-liquid fixed point theory with the limit A/kp — > 0. In 
sum, we arrive at the picture where the chiral liquid solution of the quantum 
effective chiral action gives the Fermi-liquid fixed point theory. The parameters 
of the four-Fermi interactions in the phonon channel are then identified with 
the fixed-point Landau parameters. 

There are two steps to apply such scheme to nuclear /hadron phenomenology. 
The first is to derive the in-medium effective Lagrangian directly from QCD 
or from the matter-free effective Lagrangian and the second is to solve the 
built in- medium effective Lagranian. The first is very difficult. So we usually 

3 The pion will be introduced in the Section 6.1 in terms of a non-local four- Fermi 
interaction that enters in the ground state property and gives the nucleon Landau 
mass formula in terms of BR scaling and pionic Fock term. 
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go to the second step, after building the in-medium effective Lagrangian by 
reasonable guesses. We assume that the effective Lagrangian satisfies 



where S e ^ is a Wilsonian effective action arrived at after integrating out 
high-frequency modes of the nucleon and other heavy degrees of freedom. 
This action is then given in terms of sum of terms organized in chiral order in 
the sense of effective theories at low energy. 

One way to build the chiral effective Lagrangian for nuclear matter has been 
studied by Furnstahl, Serot and Tang [11,12)0 They formulated their theory 
in terms of a chiral Lagrangian constructed by using the terms which are 
governed by QCD symmetry and applying the "naturalness" condition for all 
relevant fields. In doing this, they introduced in the FTS1 a quarkonium field 
that is associated with the trace anomaly with its potential constrained by 
Vainshtein et a/.'s low energy theorem [40]. And Georgi's "naive dimensional 
analysis" [41][^] was used in the FTS2 instead of the trace anomaly. It was 
argued in [12] that a Lagrangian so constructed contains in principle arbitrarily 
higher-order many-body effects including those loop corrections that can be 
expressed as counter terms involving matter fields (e.g., baryons). This is 
essentially equivalent to Lynn's chiral effective action [39] that purports to 
include all orders of quantum loops in chiral expansion supplemented with 
counter terms consistent with the order to which loops are calculated. Though 
it is a little hard to define the fluctuation on its ground state, their models 
are very successful in describing the bulk properties of nuclei. 

Another way is to apply the strategy of the effective theory to the in-medium 
Lagrangian. The parameters of the effective Lagrangian are related to the 
vacuum state at a given density so depend on the density. One famous example 
is BR scaling [14]. Brown and Rho point [42] that the mean field solution of 
the chiral effective Lagrangian given by BR scaling approximates 



The detail of BR scaling is reviewed in the next section. The aim of this review 
is to cast BR scaling in a suitable form starting with a chiral Lagrangian de- 
scription of the ground state as specified above around which fluctuations in a 



4 The model in [11] is referred to as FTS1. That in [12] shall as FTS2. 

5 Since the strong interactions have two relevant scales /„• and A x and A x is much 
larger than f n , we can apply Georgi's naive dimensional analysis to the low energy 
hadron physics. 




(67) 



8S eff = 0. 



(68) 
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various flavour sectors are to be made. To do this we study phenomenologically 
successful FTS1 model here. 

3.4 A chiral effective model: FTS1 model 

Furnstahl, Tang and Serot [11] constructed an effective nonlinear chiral model 
that will be referred to as FTS1 model in this review. The FTS1 model incor- 
porates the scale anomaly of QCD in terms of a light ( "quarkonium" ) scalar 
field S and a heavy ( "gluonium" ) scalar field q and gives a good description of 
basic nuclear properties in mean field. One can also build a model appealing 
to general notions of effective field theories such as "naturalness condition" 
as in [12] This avoids the use of the scale anomaly of QCD. The FTS2 model 
is also an effective mean field theory which gives an equally satisfactory phe- 
nomenology as the FTS1. However the FTS1 is found to be more convenient 
for studying the role of the light scalar field in the scale anomaly we are inter- 
ested in since the FTS1 Lagrangian includes the scale anomaly term explicitly. 
In addition, Li, Brown, Lee, and Ko found that the FTS1 model reproduces 
quite successfully nucleon flow in heavy ion collisions [26]. 

The underlying assumption in the FTS1 is that the light scalar field transforms 
under scale transformation as 

S( a -V) = a d S{x) (69) 

with a parameter d that can be different from its canonical scale dimension, 
i.e. unity, while the scale dimension of the heavy gluonium, which is integrated 
out in the effective Lagrangian for normal nuclear matter, is taken to be unity. 
This assumption imposes that quantum fluctuations in the scalar channel be 
incorporated into an anomalous dimension d an = d — 1 7^ 0. A RG flow 
argument in Section 3.3 justifies this assumption heuristically. One further 
assumption of the FTS1 is that there is no mixing between the light scalar 
S(x) and the heavy scalar q in the trace anomaly. The FTS1 Lagrangian has 
the form 

r// ^-"4( in H)-Mf) ! (>f4) <™> 

where C s is the chiral- and scale-invariant Lagrangian containing q, S, N, ir, u>, 
etc. Here ^ and S are the vacuum expectation values(VEV) with the vacuum 
|0) defined in the matter- free space: 

*> = <0k|0), S = (0\S\0). (71) 
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The trace of the improved energy- momentum tensor [37], i. e. the divergence 
of the dilatation current _D M , from the Lagrangian is; 



'S 



2\ 2 /d 



(72) 



The mass scale associated with the gluonium degree of freedom is higher than 
that of chiral symmetry, A x ~ 1 GeV. For example the mass of the scalar 
glueball calculated by Weingarten [43] is 1.6~1.8 GeV. This invites us to 
integrate out the gluonium. The resulting FTS1 effective Lagrangian takes 
the form 



£ = N[i llx {d» + iv" + igyu" + gAlsa?) -M + g s <j>]N 



+ \ (l + |f tr (dpUPUi) + m 2 v u,u^ 



m 



-a 



4,/d 



V \ So, 



11 1 

4 K 



(73) 



where S = S — 0, r/ and ( are unknown parameters to be fixed and 



£ 2 = jj = e i7r ' T I** 

v^ = ~(^ + ^) 

Note that a given VEV of the field scales down the pion decay constant and 
the uj mass in the same way at the lowest chiral order as in [14]. The static 
mean field equations of motion for FTS1 are 



i-d 

9s £ = V 2 O - ml (l - |) " In (l - |) - ^-m^ (74) 

g v £ NjN t = V 2 cu + m 2 v (l + u + ^> 3 (75) 

with the static mean field solutions O an d c^o- Equation (75) is a constraint 
because u is not a dynamical degree of freedom. 

It is important to note that the FTS1 Lagrangian is an effective Lagrangian in 
the sense explained in Section 3.2. The effect of high frequency modes of the 
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nucleoli field and other massive degrees of freedom appears in the parameters 
of the Lagrangian and in the counter terms that render the expansion mean- 
ingful. It presumably includes also vacuum fluctuations in the Dirac sea of the 
nucleons [11,13]. In general, it could be much more complicated. Indeed, one 
does not yet know how to implement this strategy in full rigor given that one 
does not know what the matching conditions are. In [11,12], the major work 
is, however, done by choosing the relevant parameters of the FTS1 Lagrangian 
to fit the empirical informations. 

The energy density for uniform nuclear matter obtained from (73) is 



kp 

e= 7 



(2nY 

m l (-1 , «M 2 , C 4 4 

— 2~ I V ~S~ j U ° 9vpBLU ° ~~ 4!^° 
Here 7 is the degeneracy factor. 



3.5 Anomalous dimension in FTS1 model 



The best fit to the properties of nuclear matter and finite nuclei is obtained 
with the parameter set Tl 6 when the scale dimension of the scalar S is near 
d = 2.7 in the FTS1 [11]. The large anomalous dimension means that one 
is fluctuating around a wrong ground state. Brown-Rho scaling is meant to 
avoid this. In this section, we analyze how this comes out and present what 
we understand on the role of the large anomalous dimension d an = rf-l«1.7 
in nuclear dynamics. In what follows, the parameter Tl with this anomalous 
dimension will be taken as a canonical parameter set. 

3.5.1 Scale anomaly 

Following Coleman and Jackiw [37], the scale anomaly can be discussed in 
terms of an anomalous Ward identity. Define T^ u (p, q) and T(p, q) by 

G( P )T llu (p,q)G(p + q) = J d 4 xd*ye^e* p -y(0 | T*9 til/ (x)(p(y)(p(0) | 0) (77) 

6 Explicitly the Tl parameters are: gl = 99.3, gl = 154.5, n = —0.496, and Q = 
0.0402. 
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d 4 xd 4 ye iq - x e ip - y (0 | T ¥ [d li D li ](x)<p(y)<p(0) | 0) 



with the renormalized propagator G(p) and the renormalized fields ip(x). Here 
T* represents the covariant T-product, D^(x) the dilatation current, and 6^ 
the improved energy momentum tensor with = x v Q iiV . A naive considera- 
tion on Ward identities would give 



with d the scale dimension of <p(x). However T is ill-defined due to singularity 
and so has to be regularized. With the regularization, the Ward identity reads 



where the additional term, A, is the anomaly. This anomaly corresponds to 
a shift in the dimension of the field involved at the lowest loop order but at 
higher orders there are vertex corrections. One obtains however a simple result 
when the beta functions vanish at zero momentum transfer [37]. Indeed in this 
case, the only effect of the anomaly will appear as an anomalous dimension. 
In general this simplification does not occur. However it can take place when 
there are nontrivial fixed points in the low energy theory. Under the reasoning 
developed in condensed matter physics [8] it is argued later that nuclear mat- 
ter is given, in the absence of BCS channel, by a Landau Fermi-liquid fixed 
point theory with vanishing beta functions of the four-Fermi interactions and 
that all quantum fluctuation effects would therefore appear in the anomalous 
dimension of the scalar field S. That nuclear matter is a Fermi-liquid fixed 
point seems to be well verified at least phenomenologically as seen in Sect. 5. 
However that fluctuations into the scalar channel can be summarized into an 
anomalous dimension is a conjecture that remains to be proved. We conjecture 
here that this is one way we can understand the success of the FTS1 model. 

3.5.2 Nuclear matter properties at d an ^5/3 

The FTS1 theory has some remarkable features associated with the large 
anomalous dimension. Particularly striking is the dependence on the anoma- 
lous dimension of the compression modulus K and many-body forces. 

In Table 1 are listed the compression modulus K and the equilibrium Fermi 
momentum k eq vs. the scale dimension d of the scalar field 0. As the d increases, 
the K drops very rapidly and stabilizes at K ~ 200 MeV consistent with 
experiments for d 2.6 and stays nearly constant for d > 2.6. This can be 



9^(p, q) = T{p, q) - idG~\p) - idG'^p + q) 



(78) 



9^(p, q) = T(p, q) - idG'^p) - idG~\p + q) + A(p, q) 
A{p,q)= lim Y{p, q, A)- Y{p, q) 



(79) 
(80) 
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Table 1 

Equilibrium Fermi momentum k eq and binding energy B = M — E/Aasa function 
of d for Fig. 4 



d K(MeV) k eq {MeV) B{MeV) 



2.3 


1960 


313 


50.4 


9 4 


1 97^ 


ouo 


^7 n 

O I .u 


2.5 


687 


297 


27.1 


2.6 


309 


279 


20.4 


2.7 


196 


257 


16.4 


2.8 


184 


241 


14.0 


2.9 


180 


231 


12.4 


3.0 


175 


223 


11.2 


3.1 


169 


217 


10.3 



seen in Fig. 4. The equilibrium Fermi momentum on the other hand slowly 
decreases uniformly as the d increases. 

Unfortunately, we have no simple understanding of the mechanism that makes 
the compression modulus K stabilize at the particular value d an ~ 5/3. We 
believe there is a non-trivial correlation between this behavior of K and the 
observation made below that the scalar logarithmic interaction brought in by 
the scale anomaly is entirely given at the saturation point by the quadratic 
term at the same d an with the higher polynomial terms (i.e. many-body inter- 
actions) contributing more repulsion for increasing anomalous dimension. At 
present our understanding is purely numerical and hence incomplete. 

In mean field, the logarithmic potential in Eqs. (73) and (76) contains n-body- 
force (for n > 2) contributions to the energy density. For the FTS1 parameters, 
these n-body terms are strongly suppressed for d > 2.6. This is shown in Fig. 
5 where it is seen that the entire potential term is accurately reproduced by 
the quadratic term \m 2 (f> 2 for d an ~ 5/3. Furthermore a close examination of 
the results reveals that each of the n-body terms are separately suppressed. 
This phenomenon is in some sense consistent with chiral symmetry [2] and is 
observed in the spectroscopy of light nuclei [44]. 

Since there are additional polynomial terms in vector fields (i.e. terms like 
(f)uj 2 ), the near complete suppression of the scalar polynomials does not mean 
the same for the total many-body forces. In fact we should not expect it. To 
explain why this is so, we plot in Fig. 6 the three-body contributions of the 3 
and 4>uj 2 forms. We also compare the FTS1 results with the FTS2 results that 
are based on the naturalness condition. In FTS1 the 3 term that turns to 
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Fig. 4. Compression modulus vs. anomalous dimension. The parameter set used 
here is the Tl in FTS1. This shows the sensitivity of the compression modulus to 
the anomalous dimension. 

repulsion from attraction for d > 8/3 contributes little, so the main repulsion 
arises from the 4>u 2 -type term. This, together with an attraction from a uo A 
term, is needed for the saturation of nuclear matter at the right density. This 
raises the question as to how one can understand the result obtained by Brown, 
Buballa, and Rho [45] where it is argued that the chiral phase transition in 
dense medium is of mean field with the bag constant given entirely by the 
quadratic term ~ \m 2 s <p 2 '. The answer to this question is as follows. First we 
expect that the anomalous dimension will stay locked at d an = d — 1 ~ 5/3 
near the phase transition (this is because the anomalous dimension associated 
with the trace anomaly - a consequence of ultraviolet regularization - is not 
expected to depend upon density), so the 0™ terms for n > 2 will continue 
to be suppressed as density approaches the critical value. Secondly near the 
chiral phase transition the gauge coupling of the vector meson will go to zero 
in accordance with the Georgi vector limit [46], where chiral symmetry is 
restored by m p —> 0, and vector meson decoupling takes place as argued in 
[47]. So the many-body forces associated with the vector mesons will also be 
suppressed as density increases to the critical density. 

3.5.3 Anomalous dimension and the scalar meson mass 

We would like to understand how the large anomalous dimension needed here 
could arise in the theory and its role in the scalar sector. As suggested in 
[21] and elaborated more in Section 3.3, one appealing way of understanding 
the FTS1 mean field theory is to consider all channels to be at Fermi-liquid 
fixed points except that because of scale anomaly, the scalar field develops an 
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Fig. 5. Comparison between the <f) interaction and the logarithmic self-interaction 
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Fig. 6. The 3-body contributions to the energy per nucleon vs. Fermi momentum in 
the FTS models. The short-dashed line represents the contribution of the <fi 3 term in 
the FTS2 with the Ql parameters. The long-dashed and the solid lines represent the 
contributions of the cubic terms (4>lo 2 and ^> 3 ) in the FTS1 with the Tl parameters 
for d = 2.7. 

anomalous dimension, thereby affecting four-Fermi interaction in the scalar 
channel resulting when the scalar field is integrated out. If the anomalous di- 
mension were sufficiently negative so that marginal terms became marginally 
relevant, then the system would become unstable as in the case of the NJL 
model or superconductivity, with the resulting spontaneous symmetry break- 
ing. However if the anomalous dimension is positive, then the resulting effect 
will instead be a screening. The positive anomalous dimension we need here 
belongs to the latter case. We can see this as follows. Consider the potential 
given with the low-lying scalar S (with the gluonium component integrated 
out): 



V( 5 ,...) = i m |^(|) 3 Q.„|-i) + ... (8!) 

where ms is the light-quarkonium mass in the free space (~ 700 MeV) and 
the ellipses stand for other contributions such as pions, quark masses etc. that 
we are not interested in. The scalar excitation on a given background S* is 
given by the double derivative of V with respect to S at S — S* 











(I) 3 " 2 


1 + 















(82) 



We may identify the ratio S* / S with the BR scaling factor $ [21]: 
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Then we have 




(84) 



m s 



with 



K d (p) = ^~ 2 l + (--l)ln$ 



(85) 



One can see that for d — 1 which would correspond to the canonical dimension 
of a scalar field the scalar mass falls much faster, for a $(p) that decreases 
as a function of density, than what would be given by BR scaling. Increasing 
the d (and hence the anomalous dimension) makes the scalar mass fall less 
rapidly. With d ~ 2, Kd ~ 1 and we recover the BR scaling. Since the drop- 
ping scalar mass is associated with an increasing attraction, we see that the 
anomalous dimension plays the role of bringing in an effective repulsion. One 
may therefore interpret this as a screening effect of the scalar attraction. In 
particular, that d — 2 pa .7 > means that in FTS1, an additional screening 
of the BR-scaled scalar exchange (or an effective repulsion) is present. 



4 Brown-Rho scaling 

Brown and Rho develop a strategy of density-dependent effective field the- 
ory for an in- medium effective theory in [14]. They assume that the effec- 
tive Lagrangian even for the hadrons in matter also keeps the symmetries of 
QCD (e.g., chiral symmetry) and that the parameters of the effective the- 
ory are determined at a given density. The change of vacuum in the presence 
of medium is assumed to be expressible by the change of parameters of the 
theory. Brown-Rho scaling is the relation among the density-dependent pa- 
rameters in medium. The quasiparticle picture associated with BR scaling is 
a successful way to describe hadrons in medium. In this section BR scaling is 
derived from QCD-oriented effective Lagrangian and discussed. 

4-1 Freund-Nambu model 

Before discussing the effective theory for dense matter and BR scaling, we 
study in this section the Freund-Nambu (FN) model, in order to describe the 
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idea of BR scaling. The FN model is a simple model where scale symmetry is 
realized in Goldstone mode. The FN model Lagrangian is 



£™ = £™ + C™ (86) 

with matter field and dilaton field x both of which have scale dimension 1; 

5<p = e(l + x ■ d)(j) 5x = e(l + x ■ d)x- (87) 

Equation (87) comes from $(ax) = a _1 $(x) with a = 1 + e where $ represents 
fields of scale dimension I. The equations of motion for the FN model 

5 2 0-/V0 = O (88) 
d 2 X - fxtf ' \xbd - v 2 ) = (89) 

have non-trivial constant solutions; 

= 0, X = ±v. (90) 
Let (x) = v and shift the dilation field as 

X = X+v. (91) 

The scale transformation of the newly defined x' does not show a definite scale 
dimension, but show a symptom of Goldstone mode; 

$ x ' = e (l +x-d)x' + ev. (92) 

With the field redefinition 



Eq. (86) becomes 
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4™ = ^ W 




1 + 2xi 'v 



1 



(94) 



with 



m ( p = f(x) 



(95) 



^ CT = v / c(x)- 
Scale symmetry in C inv becomes invisible. 

We can see some interesting features in the FN model. Both the matter and 
dilaton field masses depend on (x). It means that the both masses move in 
the same way if the vacuum expectation value of x is changed. In addition, 
the a mass can be left massive even if c — > 0. We must also note that the 
explicit scale symmetry breaking term is necessary for the FN model picture. 
It makes the spontaneously scale symmetry breaking of C inv possible. 

4-2 Brown-Rho scaling in in-medium effective Lagrangian 

We can introduce the elementary scalars that QCD has. As Schechter did [48], 
the scalar fields represented by the trace anomaly (60) can be introduced. 
We divide the trace anomaly into a hard part, which is associated with the 
gluonium Xg, an d a soft part, which comes from the quarkonium x- 



The gluon condensate in (9^)h ar d determines the gluonium mass m Xg ~ 1.6 — 
1.8 GeV [43] and does not vanish [49]. Since gluonium mass is much larger 
than the chiral scale A x ~ 1 GeV, gluonium fields are integrated out in low 
energy physics. And Adami and Brown [50] show by QCD sum rule that the 
gluon condensate is less important for the masses of light-quark hadrons. So 
we focus on the quarkonium scalar which has mostly to do with (Ofysoft- 

In matter-free space, there seems to be no scalar whose mass is small enough, 
i.e.,<^C A x . However, according to Weinberg's mended symmetry [51] there 
must be a scalar to form a quartet with pions near the chiral phase transition. 
In addition, lattice simulations [52] near the chiral phase transition shows that 
two light flavor QCD transition reproduces a scaling relation with 0(4) expo- 
nents as argued first by Pisarski and Wilczek [53]. Beane and van Kolck [54] 
suggest that the Goldstone boson from spontaneous scale symmetry breaking 



^ - (9%) hard + (0£)soft- 



(96) 
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plays the role of the chiral partner of the pion, i.e., the chiral singlet scalar 
in the scale anomaly approaches the pions and makes up the quartet of 0(4) 
symmetry in medium as density increases to the critical density of the chiral 
restoration. And we must note that the quark condensate, which contributes 
to the quarkonium mass, goes to zero as chiral symmetry is restored. So we 
assume that the quarkonim mass in the scale anomaly becomes m* <C A x as 
density increases, though m x ~ A x in free space. 

Integrating out the hard part of the scalar field Xg, we introduce into a Skyrme 
Lagrangian - which represents the low energy QCD for the infinite number 
of colors - the soft part of scalar, x, in order to make our effective theory 
consistent with QCD scale property. 



C = C inv + Csb (97) 
= f Tr (d,Ud^) + \d^X 

+ ^Tr [tfd^U, U^d v U] 2 + ■■■ (98) 
Csb = — V{x) + pion mass term + • • • (99) 

with 



\J = e V- (100) 
Xo = (0|x|0). (101) 

We make Ci nv scale correctly by multiplying the proper powers of X- Since 
Ci nv do not contribute to 6fc, we add the scale breaking potential term V(x) 
due to scale anomaly. V(x) includes radiative corrections of high order and 
gives the trace anomaly of QCD in terms of x- 



Let us break the scale symmetry spontaneously following the strategy of the 
FN model. If we define x* as the mean field value of x m dense matter 



X, = (X) P , (102) 
we can expand x as 



X = X, + X'- (103) 

Stars will be affixed to all the quantities that appear in nuclear matter from 
now on. C inv becomes in terms of x', 
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+ ^Tr [U%U, U^d v U? + ■■■ (104) 
with the effective pion decay constant defined at mean field level as 



Xo 

Note that U* = e' T-7r */^ with ir* = nx*/Xo ls the same as U in (100) by the 
definition (105). In our effective field theory both scale symmetry and chiral 
symmetry are realized in the Goldstone mode and their Goldstone bosons are 
X and 7r's, respectively. 

According to Gell-Mann-Oaks-Renner (GMOR) relation 



/X = -^^(0\uu + dd\0), (106) 

The pion mass is proportional to the quark mass. Since the quark masses comes 
from explicit chiral symmetry breaking which has something to do with the 
electroweak symmetry breaking scale Aew ^ A x , the pion mass problem is 
out of our interesting range. In the prediction of chiral perturbation theory 
in medium [55] the pole mass of the pion does not decrease up to nuclear 
matter density. In fact a recent analysis of deeply bound pionic states in heavy 
nuclei [56] shows that the pole mass of the pion could be even a few per cents 
higher than the free space value at nuclear matter density. The m* in our 
in-medium effective chiral Lagrangian is not necessarily the pole mass. Thus 
it is not clear how to incorporate this empirical information into our theory. 
We will assume here that m* does not scale. This assumption may be justified 
by using the fact that the pion is an almost perfect Goldstone boson. Under 
the assumption m* = m n , Eq. (106) may give 



Hlrf- 

It is an example to show the relation between the BR scaling factor and the 
change of the vacuum in medium. 

Since the hedgehog solution of the Skyrme Lagrangian (98) gives the Skyrmion 
mass proportional to J^f-m the nucleon mass must scale in the matter as 



42 



m* 



N 



M 




(108) 



Here M represents the nucleon mass in free space and m* N is the in-medium 
effective mass of nucleon (later identified with the Landau effective mass). In 
the simple Skyrme model g& is inversely proportional to e 2 which does not 
scale since the quartic Skyrme term is classically scale-invariant. So does 
not change in our simple model and the nucleon mass scales in the same way 
as pion decay constant in medium. 

The successful low energy results of the tree level in the chiral effective the- 
ory implemented with hidden local symmetry, i.e., Kawarabayashi-Suzuki- 
Riazuddin-Fayyazudin(KSRF) relation, p coupling universality, p-meson dom- 
inance, etc., are shown to remain valid with the loop effects at low energy [57]. 
Though it is proven in free space, it might hold in medium. If we assume KSRF 
relation 

m 2 v = 2e 2 fl (109) 
is satisfied in nuclear matter, we obtain 



m* f* 

— ^ = -j- (110) 

with the subscript v standing for light-quark vector mesons p and uj. 

The fluctuating component x' m f ne s °ft X is expected to represent multi- 
pion excitations in scalar channel and to give a scalar effective field a in dense 
medium by interpolation. It is known that the correlated 2n exchange can 
be approximated by a scalar field with a broad mass distribution [58] . Durso, 
Kim, and Wambach's recent calculation of NN — > tttt helicity amplitude in 
the scalar channel f^ =0 for p ~ p with BR scaling of vector meson mass shows 
that the resonance of the scalar mass becomes very sharper and that its value 
shifts downward 500 MeV [59]. p represents the normal nuclear matter 
density (0.16/fm 3 ). The light (<C A x ) and decreasing mass of the scalar field 
suggests that it is the expected Beane and van Kolck's dilaton [54] . Since the 
main 7r — n rescattering comes from p-meson exchange [58], the shift of the 
scalar mass in medium is affected by the density-dependence of the p mass. It 
is clear for low densities, where the linear approximation works well, that 



<^<. (Ill) 
m a m v 

Now we find that the hadron masses and pion decay constant decrease in the 
similar way. 
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Hp) 



fM 



m*(p) 



m*(p) 



m*(p) 
M 



(112) 



m, 



■v 



For the moment, we are ignoring the scaling of Qa to which we will return 
later. M* represents the effective nucleon mass obtained with a fixed g A . The 
universal factor $(p) can be determined by experiments and/or QCD sum 
rules. 

Note that the scalar field that governs BR scaling is the quarkonium compo- 
nent of the scale anomaly, not the hard gluonic component. The latter gives 
dominant contribution to the scale anomaly in QCD but is integrated out in 
the effective Lagrangian. This structure imposes the hadron scaling relation 
(112) by a Nambu-Jona-Lasinio(NJL) mechanism as described in [45]. Re- 
cently Liu et aVs detailed lattice analysis [60] for the source of the mass of 
a constituent quark supports this structure. They show that dynamical sym- 
metry breaking contribution gives most of the mass of the chiral constituent 
quark. It means that the change of the vacuum in (112) can be related only 
in a subtle way to the light quark hadron mass. 

When considering BR scaling, one must note that the BR-scaled masses and 
decay constants are the parameters in an effective theory. An effective pa- 
rameter defined in one theory can be different from the parameter defined in 
another theory, even though two theories have the same physics in common. So 
the connection between the BR-scaled parameters and the parameters in the 
other theories needs be established via observed quantities, i.e., experiments. 

One must also note that BR scaling is the result of mean field approximation. 
So when one considers the cases where higher order modifications are impor- 
tant, that is, the mean- field approximation is not reliable, BR scaling cannot 
be applied without further corrections. Let us take Goldberger-Treiman rela- 
tion 



as an example. In the real world gA decreases from ~ 1.26 in free space to 
~ 1 at normal nuclear matter density since it is affected by the short range 
interaction between baryons. So the naive BR scaling, which would have g n NN 
remain constant, can not explain Goldberger-Treiman relation in the low den- 
sity region. It means higher order modifications spoil the tree-order BR scaling 
at low densities. At high densities above the normal nuclear matter density, 
however, g\ remains at 1 while /* continues to drop and hence the coupling 
constant ratio increases. 

Although Brown and Rho's arguments about the scaling relation of effec- 
tive parameters may seem a bit drastic, the experiments (e.g., CERES [15], 



g^NNfn = gAM 



(113) 
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central S+Au, 200 GeV/nucleon _ 




M (MeV) 



Fig. 7. The comparison of CERES dilepton experiments and the theoretical predic- 
tions with the free- mass mesons and with scaled mass mesons. The figure comes 
from [61]. 

HELIOS-3 [16]) provide support for this scheme. The explanation of CERES 
data is one of the good examples. As seen in Fig. 7, the scaling of effective 
masses of hadrons in medium reproduces the CERES results very simply at 
the mean field level [17] .Q 

4-3 Duality 

Brown-Rho scaling describes the low-mass dilepton enhancement on the ba- 
sis of the quasiquark picture. But the increase of dilepton yields can also be 
described by a hadronic theory with free masses, introducing the appropri- 
ate variables, e.g. nucleonic excitations. Rapp, Chanfray and Wambach [63] 
described the CERES results successfully using conventional many-body the- 
ory. To evaluate the in-medium rho meson propagator they renormalized the 
intermediate two-pion states, which dress rho mesons and interact with the 
surrounding nucleons and deltas, and considered the direct interaction of rho 

7 In Ref. [17] BR scaling is realized on the basis of the constituent quark model [62] . 
Li, Ko, and Brown used Walecka theory to implement the density dependence of 
the constituent quarks. 
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mesons with surrounding hadrons, especially p-like baryon-hole excitations 
( "rhosobar" ) . They found that such medium effects broaden the spectral func- 
tion of rho meson and the dilepton production at ~ m p /2 is enhanced. 

Rapp et al.'s hadronic rescattering approach and Brown- Rho's quasiparticle 
approach have merits and demerits, compared with each other. Since BR scal- 
ing is approximate at mean-field level, Rapp et a/.'s theory is more reliable at 
low densities than BR scaling. On the other hand at higher densities many 
diagrams have to be considered in the hadronic rescattering approach. BR 
scaling gives a medium-modified vacuum around which the weak fluctuations 
can be dealt with. 

If both descriptions are correct, the effective variables are expected to change 
smoothly from hadrons to quasiquarks subject to BR scaling and both de- 
scriptions must show duality around the hadron-quark changeover densities^ 
Such duality was suggested by Brown et al. [65] and Y. Kim et al. [66] studied 
it more precisely. For rho mesons in medium they studied a two-level model, 
which consists of the collective rhosobar [N* (1520) A r_1 ][^| and the 'elementary' 

P- 

The in-medium p-meson propagator is 



(114) 

where the £ indicates self-energies and m° is the bare mass of p. Taking the 
free rho meson mass m p = {mP p ) 2 + Re S p7r7r , we obtain the dispersion relation 
in the q = limit 



D r 



K) 2 



J p7T7T 



— £ 



pBN 



q =m p + ReY: pN , N . (115) 

The phenomeological Lagrangian for the s-wave interaction^] between the 
elmentary p and the p-like baryon-hole excitation is 



8 Rapp and Wambach suggested /citerw99 recently interpreting the strong broad- 
ening of the /3-meson resonance as a reminiscence of hadron-quark changeover. 

9 iV*(1520) gives the most important contribution to the photoabsorption cross 
sections in the dileption analysis [67]. 

10 For p-wave coupling, the Lagrangian is 

C = &±B\s x p) ■ p a t a N + h.c. 
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Fig. 8. The in-medium p-meson mass and Z factor obtained in [66] for F t = 0. 
C p bn = fj ^B\q Q s ■ Pa - p° a s ■ q)t a N + h.c. (116) 



VClp 



with appropriate spin operator s and isospin operators t a [68] . The self-energy 
from the interaction (116) for B = A r *(1520) is 



^omsigo) ~ jW^ T ^(go + iT t /2) 2 - (AE) 2 ) (117) 

where p n is the nuclear density and the total width T t includes the the free 
width of M N * (1520) and its modification in medium. Neglecting nuclear Fermi 
motion (q = limit), AE = M N * - M N . 

Kim et al. [66] showed that the dispersion relation (115) gives the two states, 
which have p-meson quantum numbers, with the spectral weight 



Z= (l-^-ReE^r^j (118) 

and that one of them can be interpreted as an in-medium vector meson whose 
mass decreases. Figure 8 shows the results with T t = 0. R/W indicates that m p 
in (117) is the free mass and B/R indicates that qo/m p in (117) is replaced by 1, 
i.e., replacing m p by to*. Note that go i-e., the in-medium p-meson mass, cannot 
go to zero at any density as seen in R/W of Fig. 8 if m p in the denominator 
of (117) is fixed. It matches neither BR scaling [14] nor the prediction that 
to* — > at the chiral transition point [45]. Brown et al. [65] suggested the 
replacement of m p in (117) by to* in order to go from Rapp et aZ.'s hadronic 
theory to BR scaling which predicts zero vector meson mass at some high 
density. B/R of Fig. 8 shows that in- medium p-meson mass goes to zero near 
p ~ 2.75po as predicted in [45]. The study of the schematic model [65,66] 
provides how BR scaling can be interpolated from a hadronic rescattering 
description. 
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5 Effective Lagrangian with BR scaling 

If the large anomalous dimension of the scalar field in FTS1 is a symptom 
of a strong-coupling regime, it suggests that one should redefine the vacuum 
in such a way that the fluctuations around the new vacuum become weak- 
coupling. This is the basis of the BR scaling introduced in [14]. The basic 
idea is to fluctuate around the "vacuum" defined at p « p characterized 
by the quark condensate (qq) p = (qq)*- This theory was developed with a 
chiral Lagrangian implemented with the trace anomaly of QCD as seen in the 
last section. The Lagrangian used was the one valid in the large iV c -limit of 
QCD and hence given entirely in terms of boson fields from which baryons 
arise as solitons (skyrmions): Baryon properties are therefore dictated by the 
structure of the bosonic Lagrangian, thereby leading to a sort of universal 
scaling between mesons and baryons. One can see, using a dilated chiral quark 
model, that BR scaling is a generic feature also at high temperature in the 
large N c limit [69] Q 

In this description, one is approximating the complicated strong interaction 
process at a given nuclear matter density in terms of "quasiparticle" excita- 
tions for both baryons and bosons in medium. This means that the properties 
of fermions and bosons in medium at p ~ po are given in terms of tree diagrams 
with the properties defined in terms of the masses and coupling constants uni- 
versally determined by the quark condensates at that density. 

The question then is: How can one marry the FTS1 Lagrangian with the BR- 
scaled Lagrangian? The next question is how to identify BR-scaled parameters 
with the Landau parameters. In this section we will provide some answers to 
these two questions. 



5.1 A Hybrid BR-scaled model 

As a first attempt to answer these questions, we consider the hybrid model 
in which the ground state is given by the mean field of the FTS1 Lagrangian 
Cftsi and the fluctuation around the ground state is described by the tree 
diagrams of the BR-scaled Lagrangian AC, 

C eff = C FT si + AC. (119) 



11 If one calculate f£(T) and m* N (T) for zero density in chiral perturbation theory, 
the temperature-dependence deviates from BR scaling at low T [70]. In [69] Y. Kim, 
H.K. Lee and M. Rho discussed the point that BR scaling holds at non-zero density. 
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This model with the canonical parameters (Tl) for the ground state and a 
BR-scaled fluctuation Lagrangian in the non-strange flavor sector was recently 
used by Li, Brown, Lee, and Ko [26] for describing simultaneously nucleon flow 
and dilepton production in heavy ion collisions. The nucleon flow is sensitive 
to the parameters of the baryon sector, in particular, the repulsive nucleon 
vector potential at high density whereas the dilepton production probes the 
parameters of the meson sector. With a suitable momentum dependence im- 
plemented to the FTS1 mean field equation of state, the nucleon flow comes 
out in good agreement with experiments. Furthermore the scaling of the nu- 
cleon mass in the FTS1 theory in dense medium, say, at p ~ 3p , is found 
to be essentially the same as that given by the NJL model. Therefore we can 
conclude that the nucleon in FTS1 scales in the same way as BR scaling. 

The dilepton production involves both baryon and meson properties, the for- 
mer in the scaling of the nucleon mass and the latter in the scaling of the 
vector meson (p) mass. The equation of state correctly describing the nucleon 
flow and the BR-scaled vector meson mass is found to fit the dilepton data 
equally well, comparable to the fit obtained in [17] using Walecka mean field. 
What is important in this process is the scalar mean field which governs the 
BR scaling and hence the production rate comes out essentially the same for 
FTS1 and Walecka mean fields. The delicate interplay between the attraction 
and the repulsion that figures importantly in the compression modulus does 
not play an important role in the dilepton process. 

Let us see how the particles behave in the background of the FTS1 ground 
state given by Cftsi- The nucleon of course scales a la BR as mentioned 
above. We can say nothing on the pion and the p meson with the FTS1 theory. 
However there is nothing which would preclude the p scaling a la BR and the 
pion non-scaling within the scheme. What is encoded in the FTS1 theory is 
the behavior of the u and the scalar S which figure importantly in Walecka 
mean fields. Let us therefore focus on these two fields in medium near normal 
nuclear matter density. 

We have already shown in subsection 3.5.3 that the mass of the scalar field 
S drops less rapidly than BR scaling for d > 2. One can think of this as a 
screening of the four-Fermi interaction in the scalar channel that arises when 
the scalar meson with the BR-scaled mass is integrated out. This feature 
and the property of the iv field can be seen by the toy model of the FTS1 
Lagrangian (that includes terms corresponding up to three-body forces) 



UoyFTSi - L-br + —{2 + r))—u — — {120) 

where 
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C B r = N{i ltl {d» + ig v u») -M + g s 6)N 

m 2 , 26. m 2 Q ,,. 26 . . 

+ t^ (1 " v ~ "f ( ~ si; } - ( } 

We have written C BR such that the BR scaling is incorporated at mean field 
level as 



*(A0 = ^ = ^ = ^«l-# (122) 
M m s m w So V ; 



with 



5 = (0\S\0) = M/g s . (123) 

Here we are ignoring the deviation of the scaling of the effective nucleon mass 
(denoted later as m* N ) from the universal scaling $(p) [21]. This will be incor- 
porated in the next subsection. We can see from (120) that the FTS1 theory 
brings in an additional repulsive three-body force coming from a cubic scalar 
field term while if one takes rj = —2, the u field will have a BR scaling mass 
in nuclear matter. Fit to experiments favors r\ —1/2 instead of —2, thus 
indicating that the FTS1 theory requires a many-body suppression of the re- 
pulsion due to the u exchange two-body force. (In the simple model with BR 
scaling that we will construct below, we shall use this feature by introducing a 
"running" vector coupling g* that drops as a function of density.) The effective 
uj mass may be written as 



< 2 ~[l+^rK- (124) 

For r] < 0, we have a falling u mass corresponding to BR scaling (modulo, of 
course, the numerical value of 77). In FTS1, there is a quartic term ~ u 4 which 
is attractive and hence increases the u mass. In fact, because of the attractive 
quartic u term, we have 



^ 1.12 (125) 

m w 

at the saturation density with Tl parameter set. This would seem to suggest 
that due to higher polynomial (many-body) effects, the lo mass does not fol- 
low BR scaling in medium. Furthermore the ui effective mass increases slowly 
around this equilibrium value. On the other hand Klingl, Waas, and Weise's 
recent sum rule analysis on current correlation function [71,72] follows BR 
scaling. It shows that effective u meson mass scales down by about 15% at 
normal nuclear matter density from its free value. We will discuss the shift of 
vector meson mass in medium in detail in Section 5.6. 



50 



5.2 Model with BR scaling 



The above hybrid model suggests how to construct an effective Lagrangian 
model that implements BR scaling and contains the same physics as the FTS1 
theory The crucial point is that such a Lagrangian is to give in the mean field 
the chiral liquid soliton solution. This can be done by making the following 
replacements in (121): 



M-g s( p ^M\ 

2/i 200 x *2 

^(l-^)-mf (126) 
and write 



C B r = N(i^(d» + ig v ^) — M* + h<t>)N 

-\n»+\^) 2 +^ 2 -^ 2 (127) 



with 



Ml = <=<= Hp) . (128) 
M m u , m 



s 



The additional term NhipN is put in to account for the difference between the 
Landau mass m* L to be given later and the BR scaling mass M*. In the chiral 
Lagrangian approach with BR scaling, the difference comes through the Fock 
term involving non-local pion exchange. This will be discussed further in the 
next chapter. For simplicity we will take the scaling in the form 



Hp) = -1 t- ( 129 ) 

l+yp/po 

with y = 0.28 so as to give $(po) — 0.78 (corresponding to kp = 260 MeV) 
found in QCD sum rule calculations [73] , as well as from the in-medium GMOR 
relation [47]. 

Note that the Lagrangian (127) treated at the mean field level would give a 
Walecka-type model with the meson masses replaced by the BR scaling mass. 
Now in order to describe nuclear matter in the spirit of the FTS1 theory, we in- 
troduce terms cubic and higher in u and fields to be treated as perturbations 
around the BR background as 
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£n-body = a(f>u 2 + b(f) 3 + au 4 + d(j) 4 + e<f) 2 uj 2 + ■ 



(130) 



where a - e are "natural" (possibly density-dependent) constants to be de- 
termined. By inserting for the and uj fields the solutions of the static mean 
field equations given by Cbr, 

mf0 = hJ2N l N l (131) 

i 

m*Jcuo = g v J2 N l N i (132) 

i 

we see that at the mean-field level, C n ^i, dy generates three- and higher-body 
forces with the exchanged masses density-dependent d la BR. Note that at this 
point, the scaling factor $ and the mean field value (131) are not necessarily 
locked to each other by self-consistency. 

As the first trial, we will consider the drastically simplified model by dropping 
the n-body term (130) and minimally modifying the BR Lagrangian (127). 
We shall do this by letting as mentioned above the vector coupling "run" 
as a function of density. For this, we use the observation made in [26] that 
the nucleon flow probing higher density requires that g^/m^ be independent 
of density at low densities and decrease slightly at high densities. We shall 
therefore take, to simulate this particular many-body correlation effect, the 
vector coupling to scale as 



q* 1 

- = ~, T- (133) 

g v 1 + zp/po 

with z equal to or slightly greater than y. The h is assumed not to scale 
although it is easy to take into account the density dependence if necessary. 



The scaling (133) seems at odds with the prediction made with the Skyrme 
model [74] where using the Skyrme model with the quartic Skyrme term in- 
versely proportional to the coupling e 2 , it was found that 




It is tempting to identify (via SU(6) symmetry) e with g v that we are dis- 
cussing here since the Skyrme quartic term can formally be obtained from a 
hidden gauge-symmetric Lagrangian by integrating out the p meson field. If 
this were correct, one would predict that the vector coupling increases - and 
not decreases - as density increases since we know that g* A is quenched in dense 
matter. This identification could be too naive and incomplete in two respects, 
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however. First of all, this skyrmion formula is a large- N c relation and secondly 
the Skyrme quartic term embodies all short-distance physics in one dimension- 
four term in a derivative expansion. Thus the constant 1/e must represent a 
lot more than just the vector-meson (p) degree of freedom. Furthermore we are 
concerned with the ui degree of freedom which in a naive derivative expansion 
would give a six- derivative term. The BR-scaled model we are constructing 
should involve not only short- distance physics presumably represented by the 
1/e term (consistent with the understanding that the quenching of is a 
short- distance phenomena) but also longer-range correlations. Therefore the 
qualitative difference should surprise no one. 

The truncated Lagrangian that we shall consider then is 
C B r = N(tl^ + ig* v u») — M* + hcf>)N 

As suggested in [75,42,5], chiral in-medium Lagrangians can be brought to 
a form equivalent to a Walecka-type model. The scalar field appearing here 
transforms as a singlet, not as the fourth component of 0(4) of the linear sigma 
model. As it stands, the Lagrangian (134) does not look chirally invariant. This 
is because we have dropped pion fields which play no role in the ground state 
of nuclear matter. In considering fluctuations around the ground state, they 
(and other pseudo-Goldstone fields such as kaons) should be reinstated. The 
chiral singlet u field and field can be considered as auxiliary fields brought 
in from a Lagrangian consisting of multi-Fermion field operators [5] via a 
Hubbard-Stratonovich transformation. 

Since we treat density-dependent parameters, we must be careful in thermo- 
dynamic consistency. After showing the way to treat density-dependent pa- 
rameters in the next section, we display the results of our model for nuclear 
matter properties. We will argue in the next subsection that the energy den- 
sity from (151) is independent of the way how the parameters move as density 
increases. 

5.3 Thermodynamic consistency in medium 

In this subsection we address the issue of thermodynamic consistency of the 
Lagrangian (134) treated in the mean field approximation. For instance, it 
is not obvious that the presence of the density-dependent parameters in the 
Lagrangian does not spoil the self-consistency of the model, in particular, 
energy-momentum conservation in the medium and also certain relations of 
Fermi-liquid structure of the matter. The purpose of this section is to show 
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that there is no inconsistency in doing a field theory with BR scaling masses 
and other parameters. This point has not been fully appreciated by workers 
in the field. The Euler-Lagrange equations of motion are in the same as the 
ones that arise from the field for the Lagrangian wherein the masses and con- 
stants are not BR scaling. While this procedure gives correct energy density, 
pressure and compression modulus, the energy-momentum conservation is not 
automatically assured. In fact, if one were to compute the pressure from the 
energy-density S, one would find that it does not give | < T u > (where T^ v 
is the conserved energy- momentum tensor and the bra-ket means the quan- 
tity evaluated in the mean-field approximation as defined before) unless one 
drops certain terms without justification. This suggests that it is incorrect to 
take the masses and coupling constants independent of fields in deriving, by 
Noether theorem, the energy-momentum tensor. So the question is: how do 
we treat the field dependence of the BR scaling masses and constants? 

One possible solution to this problem is as follows. In Section 4.2, the density 
dependence of the Lagrangian arose as the vacuum expectation value of the 
scalar field x that figures in the QCD trace anomaly. By vacuum we mean 
the state of baryon number zero modified from that of true vacuum by the 
strong influence of the baryons in the system. See later for more on this point. 
It corresponded to the condensate of a quarkonium component of the scalar 
X with the gluonium component - which lies higher than the chiral scale 
- integrated out. It was assumed to scale in dense medium in a Skyrmion- 
type Lagrangian subject to chiral symmetry. Now in the language of a chiral 
Lagrangian consisting of the nucleonic matter field N with other massive fields 
integrated out, this scalar condensate would be some function of the vacuum 
expectation value of NN or N^qN coming from multi-Fermion field operators 
mentioned above. How these four- Fermi and higher- Fermi field terms can lead 
to BR scaling in the framework of chiral perturbation theory was discussed in 
[5]. We shall follow this strategy in this paper leaving other possibilities (such 
as dependence on the mean fields of the massive mesons) for later investigation. 
For this, it is convenient to define 



(135) 



with unit fluid 4-velocity 




1 



1 



(136) 



with the baryon current density 



j =< N~fN > 



(137) 
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and the baryon number density 



p=< AftjV >=Y^ni. (138) 

i 

We will take rtj to be given by the Fermi distribution function, rij = 9(kp — \ki\) 
at T = 0. We should replace p in (134) by p for consistency of the model. The 
definition of p makes our Lagrangian Lorentz-invariant which will later turn 
out to be useful in deriving relativistic Landau formulas. With this, the Euler- 
Lagrange equation of motion for the nucleon field is 



5C dC dC dp 

+ 



5N dN dpdN 

= [«7"(fy + igfa - iujb) - M* + h<f>]N 

= (139) 



with 



dC 

dp 



(140) 



2 9r K *,2 9m t at m. *r d 9v *T*T dM 



■k 



dp s dp dp dp 



This additional term which may be related to what is referred to in many- 
body theory as "rearrangement terms" plays a crucial role in what follows. f 17 ] 
The equations of motion for the bosonic fields are 

(d^df, + mf)4> = hNN (141) 
d v F^ + mfuf = g*N^N. (142) 

We start with the conserved canonical energy-momentum tensor constructed 
a la Noether from the Lagrangian (134): 



= iNYd u N + d»<pd v <p - d»co x d u u x 

-\[(d(p) 2 - mfcf) 2 - {duf + m*V - 2±NylN]g^. (143) 

We shall compute thermodynamic quantities from (143) using the mean field 
approximation which amounts to taking 



12 For a recent discussion on rearrangement terms as well as thermodynamic consis- 
tency in the context of standard many-body theory, see [76] . The notion of density- 
dependent parameters in many-body problems is of course an old one [77]. 
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1 y^ a -J EKi + m * N 



w^r^i 2E, 

h<j> = C 2 h <NN >=C 2 h J2 



X 



\ 



(TKj 



ru 



m* N 



rn 



exp (iKi ■ x - i(g*ujQ - m S + E^t) 

(144) 



Hi 



where a$ is the annihilation operator of the nucleon i, with 



Hi =< a\ai >, 



and 



E = (E>, 

K i = ki — C 2 j + wS, 



and 



E K 



k 2 + rri$. 



X is the spinor and <x is the Pauli matrix. We have defined 



CM = 



m*(p) 



and 



(145) 



(146) 
(147) 



(148) 



(149) 



C h {p) 



h 



™*{p) C h (p) 
In this approximation, the energy density is 



(150) 



£ = <T m > 
= < iN-f°d°N + \mf<p 2 - \m**u 2 + tN^/N > 

1 ^ 



(151) 



Note that the El-dependent terms cancel out in the comoving frame (v — 0), 
so that the resulting energy-density is identical to what one would obtain from 
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the Lagrangian in the mean field with the density-dependent parameters taken 
as field-independent quantities. 

Given the energy density (151), the pressure can be calculated by (at T = 0) 



dE 2 dS/p 
= \C 2 MP 2 - S P - l& h (p)(m* N - M\p)f 

~h [ EF ^r kF - Y2 k ^ - ifM(** + ( 152 ) 

where p is the chemical potential - the first derivative of the energy density 
with respect to p in the comoving frame (v — 0): 



p = ^£\v=o = C 2 vP + Ep — S (153) 



with Ep — ykp + m*^ and S = (E) v=0 . To check that this is consistent, we 
calculate the pressure from the energy- momentum tensor (143) in the mean 
field at T = 0: 



Pt = 3 < Ta > v=0 (154) 
= ^(tN^N - \{m^u? - mftf - 2tN^N) 9ii ) v=0 
= \C 2 MP 2 - \c 2 h {p){m* N - M*{p)f - S p 

This agrees with (152). Thus our equation of state conserves energy and mo- 
mentum. 



We showed that a simple effective chiral Lagrangian with BR scaling parame- 
ters is thermodynamically consistent, a point which is important for studying 
nuclear matter under extreme conditions. It is clear however that this does 
not require that the masses appearing in the Lagrangian scale according to 
BR scaling only. What is shown in this subsection is that masses and coupling 
constants could depend on density without getting into inconsistency with 
general constraints of chiral Lagrangian field theory. This point is important 
for applying (134) to the density regime p ~ 3po appropriate for the CERES 
dilepton experiments and also kaon production at GSI(Gesellschaft fur Schw- 
erionenforschung) where deviation from the simple BR scaling of [14] might 
occur. 
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Table 2 

Parameters for the Lagrangian (134) with y = 0.28, m s = 700MeV, m w = 783MeV, 
M = 939MeV 



SET 


h 


9v 


z 


SI 


6.62 


15.8 


0.28 


S2 


5.62 


15.3 


0.30 


S3 


5.30 


15.2 


0.31 



Table 3 

Nuclear matter properties predicted with the parameters of Table 2. The effective 
nucleon mass (later identified with the Landau mass) is m* N = M* — hfo. 



SET 


E/A — M(MeV) 


k eq (MeV) 


K(MeV) 


m* N /M 


$(k eq ) 


SI 


-16.0 


257.3 


296 


0.619 


0.79 


S2 


-16.2 


256.9 


263 


0.666 


0.79 


S3 


-16.1 


258.2 


259 


0.675 


0.78 



5.4 Results 

Based on the thermodynamic consistency of density-dependent effective the- 
ories proven at the mean field level, we check whether the model Lagrangian 
(134) can describe the infinite nuclear matter properties successfully. The char- 
acteristic properties we try to reproduce are compression modulus, m* N , and 
binding energy at normal nuclear matter density, and the saturation density 
itself. 

In Table 2, three sets of parameters are listed. We take the measured free- 
space masses for the u and the nucleon and for the scalar <p for which the 
free-space mass cannot be precisely given, we take m s = 700 MeV (consistent 
with what is argued in [47]) so that at nuclear matter density, it comes close to 
what enters in the FTS1. The resulting fits to the properties of nuclear matter 
are given in Table 3 for the parameters given in Table 2. These results are 
encouraging. Considering the simplicity of the model, the model - in particular 
with the S2 and S3 set - is remarkably close in nuclear matter to the full 
FTS1. The compression modulus comes down toward the low value that is 
currently favored. In fact, the somewhat higher value obtained here can be 
easily brought down to about 200 MeV without modifying other quantities if 
one admits a small admixture of the residual many-body terms (130), as we 
shall shortly show. The effective nucleon Landau mass m* N /M « 0.67 is in 
good agreement with what was obtained in QCD sum rule calculations [78] 
and also in the next section (i.e., 0.69) by mapping BR scaling to Landau- 
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Migdal Fermi-liquid theory. We shall see below that this has strong support 
from low-energy nuclear properties. What is also noteworthy is that the ratio 
7Z = {g^/rnlj) 2 forced upon us - though not predicted - is independent of 
the density (set SI) or slightly decreasing with density (sets S2 and S3), as 
required in the nucleon flow data as found by Li, Brown, Lee and Ko [26]. 
In FTS1 theory, it is the higher polynomial terms in uj and defining the 
mean fields that are responsible for the reduction in 1Z needed in [26]. In 
Dirac-Brueckner-Hartree-Fock theory, it is found [79] that while 1Z takes the 
free-space value 1Z for p ps* p , it decreases to 1Z ps O.QATZq at p ps 3p due to 
rescattering terms which in our language would correspond to the many-body 
correlations. 

The assumption that the many-body correlation terms in (130) can be entirely 
subsumed in the dropping vector coupling may seem too drastic. Let us see 
what small residual three-body and four-body terms in (130) as many-body 
correlations (over and above what is included in the running vector coupling 
constant) can do to nuclear matter properties. For convenience we rewrite 
(130) by inserting dimensional factors as 



r _ Vo 2 2 

1~n-body — -^ m ujJ~ UJ 



«3 2 



(155) 



Co 2 4 K 4 2 4 . Vl 2 <t> 2 2 



4! 



and demand that the coefficients rj, C and k so defined be natural. The results 
of this analysis are given in Table 4 and Fig. 9 for various values of the residual 
many-body terms and compared with those of the truncated model (134) with 
S3 parameter set. The coefficients are chosen somewhat arbitrarily to bring 
our points home, with no attempt made for a systematic fit. (It would be 
easy to fine-tune the parameters to make the model as close as one wishes 
to FTS1 theory.) It should be noted that while the equilibrium density or 
Fermi momentum k eq , the effective nucleon mass m* N and the binding energy 
B stay more or less unchanged, within the range of the parameters chosen, 
from what is given by the BR-scaled model (134) with the S3 parameters, the 
compression modulus K can be substantially decreased by the residual many- 
body terms. Figure 9 shows that as expected, lowering of the compression 
modulus is accompanied by softening of the equation of state at p > p . While 
the equilibrium property other than the compression modulus is insensitive to 
the many-body correlation terms, the equation of state at larger density can 
be quite sensitive to them. This is because for the generic parameters chosen, 
the m* N can vanish at a given density above po a t which the approximation is 
expected to break down and hence the resulting value cannot be trusted. The 
B2 and B4 models do show this instability at p > 1.5p - (See Fig. A.l.) 
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Table 4 

Effect of many-body correlations on nuclear matter properties using the La- 
grangian (134) + (155). We have fixed the free-space masses m s = 700MeV, m w = 
783MeV,M = 939MeV and set r\\ = for simplicity. The equilibrium density k eq , 
the compression modulus K, and the binding energy B = M — E/A are all given in 
units of MeV. 



SET 


h 


9v 


y 


z 


^0 


Co 


^3 


K4 


keq 


m N 

M 


K 


B 


S3 


5.30 


15.2 


0.28 


0.31 










258.2 


0.675 


259 


16.1 


Bl 


5.7 


15.3 


0.28 


0.30 






0.5 


-4.9 


256.0 


0.666 


209 


16.2 


B2 


5.7 


15.3 


0.28 


0.30 


-0.055 


0.18 






257.3 


0.661 


201 


16.1 


B3 


5.6 


15.27 


0.28 


0.30 






0.31 


-4.1 


259.1 


0.659 


185 


16.1 


B4 


5.6 


15.3 


0.28 


0.31 






0.9 


-8.1 


256.4 


0.669 


191 


16.1 


CI 


5.7 


15.3 


0.28 


0.30 


-0.05 


0.155 






256.3 


0.665 


218 


16.2 


C2 


5.8 


15.3 


0.28 


0.30 


-0.11 


0.35 






256.1 


0.662 


161 


16.2 




Fig. 9. E/A-M vs. p for FTS1 theory ("Tl" parameter), the "S3", "Bl" and "B3" 
models defined in Table 4. 
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It is quite encouraging that the simple minimal model (134) with BR scaling 
captures so much of the physics of nuclear matter. Of course, by itself, there is 
no big deal in what is obtained by the truncated model: It is not a prediction. 
What is not trivial, however, is that once we have a Lagrangian of the form 
(134) which defines the mean fields, then we are able to control with some con- 
fidence the background around which we can calculate fluctuations, which was 
the principal objective we set at the beginning of this approach. The power of 
the simple Lagrangian is that we can now treat fluctuations at higher densities 
as one encounters in heavy-ion collisions, not just at an equilibrium point. The 
description of such fluctuations does not suffer from the sensitivity with which 
the equation of state depends at p > p on the many-body correlation terms 
(130). 

The simple Lagrangian (134) embodies the effective field theory of QCD dis- 
cussed by Furnstahl, Serot, and Tang [11] anchored on general considera- 
tions of chiral symmetry. This Lagrangian should be viewed as an effective 
Lagrangian that results from two successive renormalization group "decima- 
tions" , one leading to a chiral liquid structure [39] at the chiral symmetry scale 
and the other with respect to the Fermi surface [8,9]. The advantage of (134) 
is that it can, on the one hand, be connected to Landau Fermi-liquid fixed 
point theory of nuclear matter and, on the other hand, be extrapolated to the 
regime of hadronic matter produced under extreme conditions as encountered 
in relativistic heavy ion processes. It would, for instance, allow one, starting 
from the ground state of nuclear matter, to treat on the same footing the 
dilepton processes observed in CERES experiments as explained in [17] and 
kaon production at SlS(Schwerionen-Synchrotron) energy and kaon condensa- 
tion in dense matter relevant to the formation of compact stars as discussed 
in [80]. 



5.5 Landau Fermi-liquid properties of the BR-scaled model 



The next issue we address is the connection between the mean-field theory 
of the chiral Lagrangian (5.2) and Landau's Fermi-liquid fixed point theory 
as formulated in Section 3.3. As far as we know, this connection is the only 
means available to implement chiral symmetry of QCD in dense matter based 
on effective field theory. For this, we shall follow closely Matsui's analysis of 
Walecka model in mean field [81] exploiting the similarity of our model to the 
latter. 
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5.5.1 Quasiparticle interactions 



The quasiparticle energy and quasiparticle interaction are, respectively, 
given by first and second derivatives with respect to nf 



(156) 

= (157) 



A straightforward calculation gives 



and 



+C4 ™-„ - M«)^ - Oi( roSr - M>)— - E M • ^ (158) 



dei 



°« + 4L/ *>P^T~ + fll h ^ il"a7i + °» 



<9p -Ej <9n,- dp dp 2 

+(m* N - M*U^Y + C h ^\ + 2(7^(< - M*)-^-{m* N - M*) 

-2C„-^(m« N - M*)_ - Cj-g^mJ, - M ) - d(m% - M«)_ 



v2 S , h dj 



«*-?t tfc (159) 



with Ei = yjkt + mtf. Note that C v , C h , and M* are functions of 



(p) = u p-u- j (160) 

in the mean field approximation. In arriving at (159), we have used the obser- 
vation that in the limit j — > 0, we have 



drii 

d 2 u 1 dj dj 



dnidnj p 2 drii drij ' 
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du 1 dj 

drii p drii ' 

<9 2 (p) 1 dj dj 



driidrij p drii drij 

and that if / is taken to be a function of the expectation value of p, then as 
j — > 0, we have 



9/ 9/ 9(p> 

(161) 



9nj d(p) dn 
^ dp 



d 2 f _ d 2 f d(p) d(p) + df d 2 (p) 



driidrij d(p} 2 drii drij d{p) driidrij 
_ 19/ 5j_ 
<9p 2 p <9p <9nj 9nj 



(162) 



In the absence of the baryon current, j = 0, the quantities and 
simplify to 



a * dM * - r>r, 9Ch v r,, rn x - r 2 " 1 ^ 

dm* N _ ~df ^h-Q^l^ini-^ - ^h-E~ 
drii i „ k 



1 + 

and 



(163) 



dj_ 

drii , y, n -,^ \k 



Writing in the standard way 



(164) 



r HO k ■ h 

fi = (2/ + 1) / -^Pi{^)km = M = k F ), (165) 



we see that the last term in (159) contributes to f\ and the sum of the rest at 
the Fermi surface (i.e. \kj\ = &f) to /q. So 



_ _ XkpEp 3Ep 
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+ (m^ - M*) 2 jfe 2 + + 2<7 fc ^Kr - M*)^-(m* N - M* 



dp dp 2 J dp dn 



n * * »^<9M* ~ 2 9M* 9 . . ~ 2/ . . d 2 M* 



and 



_ XkpEp 



3(g, 2 ~ g)p 



(167) 



5.5.2 Some relations for relativistic Fermi-liquid 

Here we bridge the model (134) to relativistic Fermi-liquid theory. For it we will 
show that the thermodynamic properties of any model like (134), which has 
density-dependent parameters and is Walecka-type, can be described in terms 
of relativistic Landau parameters derived from the mean field approximation 
of the model in the same way as in Section 2.3. 

First let us calculate the compression modulus K defined by 



It comes out to be 



K = S + 9p[C 2 + AC v p 8 -^ + H±^k + p 2 {(^) 2 + Cv ^} (169) 
hjp op hp op op Ofr 

Fir 1 Ffic 1 Fir 1 fi 

+Kv - M*) 2 {(^) 2 + C h d -0} + 2C h ^(m* N - M*)£ p (m* N - M*) 

-2C h —(m N — M )~q^ C ^Q-p {mN ~ M) ~ Ch{mN ~ M 
Comparing (166) and (169), we obtain Eq. (29); 



K = ^(l + F ). (170) 

Hit? 
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In our model 




E F EE Jkj 




(171) 



It is verified that our model satisfies the relativistic Landau Fermi-liquid for- 
mula for the compression modulus (29). 

As shown in Section 2.3, the relativistic Landau liquid satisfies the mass rela- 
tion (40); 



\ 1 J k=kp,v=0 

One can see from Eqs. (153), (167), and (171) that (40) is satisfied exactly in 
our model. 

And lastly the relativistic relation for the first sound velocity (46) is satisfied 
automatically since (29) and (40) are satisfied in our model. So all the rela- 
tivistic Landau Fermi-liquid relations in Section 2.3 are satisfied with density- 
dependent Walecka-type models like (134). 

5.5.3 Discussions 

The crucial question is really how to understand the scaling masses and con- 
stants as one varies temperature and density as considered in [14]. If one takes 
the basic assumption that the chiral Lagrangian in the mean field with BR 
scaling parameters corresponds to Landau's Fermi-liquid fixed point theory, 
then one should consider first fixing the Fermi momentum hp and let renor- 
malization group flow come to the fixed points of the effective mass m* L for 
the nucleon and Landau parameters T [8]. In this case, the scaling quantities 
would seem to be dependent upon A/hp, not on the fields entering into the 
effective Lagrangian. This paper however shows that if one wants to approach 
the Fermi-liquid fixed point theory starting from an effective chiral Lagrangian 
of QCD, it is necessary to take into account the fact that the scaling arises 
from the effect of multi-Fermi interactions figuring in chiral Lagrangians as 
implied by chiral perturbation theory described in [5]. This is probably due 
to the fact that we are dealing with two-stage "decimations" in the present 
problem - with the Fermi surface formed from a chiral Lagrangian as a non- 
topological soliton (i.e., "chiral liquid" [39]) - in contrast to condensed matter 
situations where one starts ab initio with the Fermi surface without worrying 
about how the Fermi surface is formed. Our result suggests that there will 
be a duality in describing processes manifesting the scaling behavior. In other 
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words, the change of "vacuum" by density exploited in [14] could equally be 
represented by a certain (possibly infinite) set of interactions among hadrons 
- e.g., four- Fermi and higher- Fermi terms in chiral Lagrangians - canonically 
taken into account in many-body theories starting from the usual matter-free 
vacuum. A notable evidence may be found in the two plausible explanations 
of the low mass enhancement in CERES dilepton yields in terms of scaling 
vector meson masses [17] and in terms of hadronic interactions giving rise to 
increased widths [63,82]. Recently the relation between two descriptions are 
discussed by Brown et al. [65] and also by Kim et al. [66] . It is argued that the 
description based on the reaction dynamics and on the meson spectral func- 
tion should be reliable at low density where the effective degrees of freedom 
are hadrons and have no contradiction with the description on BR scaling. 
There should however be a "crossover" region at higher density at which BR 
scaling will become more efficient or we should go over to constituent quarks. 
How to relate the two description in the "crossover" density regime remains 
an open problem. 

In discussing the properties of dense matter, such as the BR scaling of masses 
and coupling constants, e.g., /* , we have been using a Lagrangian which pre- 
serves Lorentz invariance. This seems to be at odds with the fact that the 
medium breaks Lorentz symmetry. One would expect for instance that the 
space and time components of a current would be characterized by different 
constants. Specifically such quantities as gu, f n etc. would be different if they 
were associated with the space component or time component of the axial cur- 
rent. So a possible question is: How is the medium- induced symmetry breaking 
accommodated in the formalism which will be discussed in the next section? 

Section 5.3 and this section provide the answer to this question. Here the 
argument is given in an exact parallel to Walecka mean field theory of nuclear 
matter. One writes an effective Lagrangian with all symmetries of QCD which 
in the mean field defines the parameters relevant to the state of matter with 
density. The parameters that become constants (masses, coupling constants 
etc.) at given density are actually functionals of chiral invariant bilinears in the 
nucleon fields. When the scalar field and the bilinear ipty, where ip is the 
nucleon field, develop a non-vanishing expectation value Lorentz invariance 
is broken and the time and space components of a nuclear current pick up 
different constants. This is how for instance the Gamow- Teller constant g& 
measured in the space component of the axial current is quenched in medium 
while the axial charge measured in the axial charge transitions is enhanced 
as described in the next chapter. If one were to calculate the pion decay 
constant in medium, one would also find that the quantity measured in the 
space component is different from the time component. The way Lorentz- 
invariant Lagrangians figure in nuclear physics is in some sense similar to 
what happens in condensed matter physics. 
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5.6 Mesons in medium 



It should be recalled that we extracted the scaling parameter $ from the in- 
medium property of the vector mesons. Here we will present evidences for the 
predicted scaling in the meson masses. There are some preliminary experimen- 
tal indications for decrease in matter of the p meson mass in recent nuclear 
experiments [83,84] but we expect more definitive results from future exper- 
iments at GSI and other laboratories. In fact, this is currently a hot issue in 
connection with the recent dilepton data coming from relativistic heavy ion 
experiments at CERN(European Organization for Nuclear Research). 

When heavy mesons such as the vector mesons p, u and the scalar a are re- 
instated in the chiral Lagrangian, then the mass parameters of those particles 
in the Lagrangian, when written in a chirally invariant way, are supposed to 
appear with star and are assumed to scale according to Eq. (112). The ques- 
tion is: What is the physical role of these mass parameters? If we assume 
that the mesons behave also like quasiparticles, that is, like weakly interact- 
ing particles with the "dropping masses," then physical observables will be 
principally dictated by the tree diagrams of those particles endowed with the 
scaling masses. In this case, the masses figuring in the Lagrangian could be 
identified in some sense as "effective" masses of the particles in the matter. 
This line of reasoning was used in the work of Li, Ko and Brown [17] to in- 
terpret the low mass enhancement of the CERES data [15]. As discussed in 
Section 5.2, this treatment is consistent with an effective Lagrangian which in 
the mean field approximation yields the nuclear matter ground state as well 
as fluctuations around the ground state. The parameters of the theory, as well 
as their density dependence are determined by the properties of the ground 
state. The work of this section shows that this scheme is internally consistent. 
However we emphasize that the scaling we have established is for the mesons 
that are highly off-shell and it may not be applied to mesons that are near 
on-shell without further corrections (e.g., widths etc.). 

Suppose one probes the propagation of an u meson in nuclear medium as in 
HADES (High Acceptance Di-Electron Spectrometer) or TJNAF (Thomas Jef- 
ferson National Accelerator Facility) experiments, say through dilepton pro- 
duction. The cu's will decay primarily outside of the nuclear medium, but let 
us suppose that experimental conditions are chosen so that the leptons from 
the ui decaying inside dense matter can be detected. See [85] for discussions 
on this issue. The question is whether the dileptons will probe the BR-scaled 
mass or the quantity given by (125). The behavior of the u mass would differ 
drastically in the two scenarios. A straightforward application of FTS1 theory 
would suggest that at a density p < po, the uj mass as "seen" by the dilep- 
tons will increase slightly instead of decrease. Since in FTS1 theory, the vector 
coupling g v does not scale, this means that (g*/m^) will effectively decrease. 
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On the other hand if the vector coupling constant drops together with the 
mass at increasing density as in the BR scaling model, the situation could be 
quite different, particularly if dileptons are produced at a density p ~ 3p as 
in the CERES experiments: The u will then be expected to BR-scale up to 
the phase transition. It has been recently suggested [86] that at some high 
density, Lorentz symmetry can be spontaneously broken giving rise to light 
u mesons as "almost Goldstone" bosons when a small explicit Lorentz sym- 
metry breaking term via chemical potential is introduced. By introducing the 
term, they assume a state which is chirally symmetric ({qq) = 0) but breaks 
Lorentz symmetry ((q^q) ^ 0). The assumed state is metastable at /i < /i c 
but becomes a global minimum at p > p c . At p > p c , u would become light 
but not massless due to the explicit breaking. Such mesons could be a source 
of copious dileptons at some density higher than normal matter density. Thus 
measuring the oo mass shift could be a key test of the BR scaling idea as 
opposed to the FTSl-type interpretations. This interesting issue will be stud- 
ied in forthcoming experiments at GSI and TJNAF. It is interesting that the 
dropping uo mass is also found in a recent QCD sum rule calculation based 
on current correlation functions by Klingl, Kaiser and Weise [71] who, how- 
ever, do not see the dropping of the p mass because of the large broadening 
of p peak. If we can describe the uj meson in medium as a quasiparticle an 
a;-nuclear bound state is feasible even in light nuclei [72]. The process like 
7 Li(d, 3 He)^,He is expected to be seen in GSI [87] if such a bound state exists. 



6 Fermi-liquid theory vs. chiral Lagrangian 

6.1 Electromagnetic current 

We will here give a brief derivation of the Landau-Migdal formula for the 
convection current for a particle of momentum k sitting on top of the Fermi 
sea responding to a slowly varying electromagnetic (EM) field. We will then 
analyze it in terms of the specific degrees of freedom that contribute to the 
current. This will be followed by a description in terms of a chiral Lagrangian 
as discussed in [21]. This procedure will provide the link between the two 
approaches. 



6.1.1 Landau-Migdal formula for the convection current 

Following Landau's original reasoning adapted by Migdal to nuclear systems, 
we start with the convection current given by 
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(172) 



where the sum goes over the spin a and isospin r which in spin- and isospin- 
saturated systems may be written as a trace over the a and r operators. More 
precisely, this is a matrix element of the current operator corresponding to 
the response to the EM field of a nucleon (proton or neutron) sitting on top 
of the Fermi sea. The sum over spin and isospin and the momentum integral 
go over all occupied states up to the valence particle. What we want is a 
current operator and it is deduced after the calculation is completed. One 
can of course work directly with the operator but the result is the same. We 
consider a variation of the distribution function from that of an equilibrium 
state 

n p = n° p + 5n p , (173) 



where the superscript refers to equilibrium. The variation of the distribution 
function induces a variation of the quasiparticle energy 



e 



p = e° p + 5e p . (174) 



In the equilibrium state the current is zero by symmetry, so we have 



J = JlJ (^3 ((^)H, + (VApK) ^(1 + rs), 

= E/ {<yv4)8np - (V P n° p )5e p )) \{l + r 3 ) (175) 

to linear order in the variation. We consider a nucleon added at the Fermi 
surface of a system in its ground state. Then 



and 



(176) 



V y p = -P-5(p-k F ) (177) 

where k with \k\ = kp is the momentum of the quasiparticle. The modification 
of the quasiparticle energies due to the additional particle is given by 
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5£ P = E / -^3fp,P ,5n p'<r'T'- (178) 

Combining (2), (175), (176) and (178) one finds that the first term of (175) 
gives the operator 



jw = Alip, (179) 

m L 2 

where k is taken to be at the Fermi surface. The second term yields 



5 J = 6J S + 5J V = h- ( El±Itl\ j (iso) 



M V 6 



where 



8J. = AIt> (181) 
m L 2 3 

S j^Jl^EL = ±^ + ±h^Ii. (182) 

ml 2 3 m£ 2 3 m\ 2 3 V ; 
For convenience, let's define^ 3 ] 

F, = (M/m* L )Fi (183) 
with analogous definitions of F/, etc.. It gives another representation of (41) 



4 = 1 - F i ■ («4) 
m L 3 

Putting everything together we recover the well known result of Migdal [31,33] 



where 



13 



The definition of m£ is in (7) 
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is the orbital gyromagnetic ratio and 

Sgi = I(^-Fi)r 3 . 



(187) 



Thus, the renormalization of Q\ is purely isovector. This is due to Galilean 
invariance, which implies a cancellation in the isoscalar channel. 

We have derived Migdal's result using standard Fermi-liquid theory argu- 
ments. This result can also be obtained [88] by using the Ward identity, which 
follows from gauge invariance of the electro-magnetic interaction. This is of 
course physically equivalent to the above formulation. We shall now identify 
specific hadronic contributions to the current (185) in two ways: the Fermi- 
liquid theory approach and the chiral Lagrangian approach. 

6.1.2 Pionic contribution 

In Fermi-liquid theory approach, all we need to do is to compute the Landau 
parameter i*\ from the pion exchange. The one-pion-exchange contribution to 
the quasiparticle interaction is 



where q = p — p' and / = g n NN(m n /2M) w 1. In a relativistic formulation 
sketched in Appendix B, we can Fierz the one-pion exchange. Done in this 
way, the Fierzed scalar channel is canceled by a part of the vector channel and 
the remaining vector channel makes a natural contribution to the pionic piece 
of Fi. The one-pion-exchange contribution to the Landau parameter relevant 
for the convection current is 



f 



P ar, P 'a'r' 3 m 2 q 2 + m 




^(S 12 (q) + 1(3-<T.<T')) 

'77- \ ^ / 



3 — T • t' 

2 



(188) 




(189) 



where 




1 



(190) 



Note that Fi(ix) satisfies 




dZ n (p) 
kp dp 



p=k F 



(191) 
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(a) 

Fig. 10. Feynman diagrams contributing to the EM convection current in effective 
chiral Lagrangian field theory. Figure (a) is the single-particle term and (b, c) the 
next-to-leading chiral order pion-exchange current term. Figure (c) does not con- 
tribute to the convection current; it renormalizes the spin gyromagnetic ratio. 

with one-pion-exchange Fock contribution to the self-energy and includes 
the higher order contribution in m* L . Thus, from Eq. (187), the one-pion- 
exchange contribution to the gyromagnetic ratio is 



Mf_ 



hr 3 - 



(192) 



In the next subsubsection we include contributions also from other degrees of 
freedom. 

Let's obtain the convection current from a chiral Lagrangian and compare it 
with the results given above. In absence of other meson degrees of freedom, we 
can simply calculate Feynman diagrams given by a chiral Lagrangian defined 
in matter-free space. Nonperturbative effects due to the presence of heavy 
mesons introduce a subtlety that will be treated below. 

In the leading chiral order, there is the single-particle contribution Fig. 10a 
which for a particle on the Fermi surface with the momentum k is given by 



Jl-body - TT (193) 

Note that the nucleon mass appearing in (193) is the free-space mass M as 
it appears in the Lagrangian, not the effective mass m\ that enters in the 
Fermi-liquid approach, (179). To the next-to-leading order, we have two soft- 
pion terms as discussed in [5,89,90,6]. We should recall a well-known caveat 
here discussed already in [90]. If one were to blindly calculate the convection 
current coming from Fig. 10b, there would be a gauge non-invariant term 
that is present because the hole line is off-shell. Figure 10c contains also a 
gauge non-invariant term which is exactly the same as in Fig. 10b but with 
an opposite sign, so in the sum of the two graphs, the two cancel exactly so 
that only the gauge-invariant term survives. Of course we now know that the 
off-shell dependence is not physical and could be removed by field redefinition 
ah initio. To the convection current we need, only Fig. 10b contributes 



72 



J2-60* = l F £~ 2 hrs = ~(^) - F l(7 r))r 3 . (194) 

We should emphasize that the Landau parameters F 1 and F[ are entirely fixed 
by a chiral effective Lagrangian for any density. 

The sum of (193) and (194) agrees precisely with the Fermi-liquid theory 
result (185) and (189). This formula first derived in [91] in connection with 
the Landau-Migdal parameter is of course the same as the Miyazawa formula 
[92] derived nearly half a century ago. Note the remarkable simplicity in the 
derivation starting from a chiral Lagrangian. However, we should caution that 
there are some non-trivial assumptions to go with the validity of the formula. 
As we will see shortly, we will not have this luxury of simplicity when other 
degrees of freedom enter. 



6.1.3 Vector meson contributions and BR scaling 

So far we have computed only the pion contribution to g t . In nuclear physics, 
more massive degrees of freedom such as the vector mesons p and uo of mass 
700 ~ 800 MeV and the scalar meson a of mass 600 ~ 700 MeV play an 
important role. When integrated out from the chiral Lagrangian, they give 
rise to effective four-Fermion interactions: 



U = %NNf - ^(A 7 ,A) 2 " ^(N^rNf + ■■■ (195) 
where the coefficients C's can be identified with 



C 2 M = JlM 1 with M = (j>,p, oo. (196) 

For the moment, we make no distinction as to whether one is taking into ac- 
count BR scaling or not. For the Fermi-liquid approach, this is not relevant 
since the parameters are not calculated. However with chiral Lagrangians, 
we will specify the scaling which is essential. Such interaction terms are "ir- 
relevant" in the renormalization group flow sense but can make crucial con- 
tributions by becoming "marginal" in some particular kinematic situation. A 
detailed discussion of this point can be found in [8]. The effective four-Fermion 
interactions play a key role in stabilizing the Fermi-liquid state and leads to 
the fixed points for the Landau parameters. (The other fixed-point quantity, 
i.e. the effective mass, is put in by fiat to keep the density fixed.) In the two- 
nucleon systems studied in [7] , they enter into the next-to-leading order term 
of the potential, which is crucial in providing the cut-off independence found 
for cut-off masses > m n . 
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Again it suffices to compute the Landau parameters coming from the velocity- 
dependent part of heavy meson exchanges in the Fermi-liquid theory approach. 
We treat the effective four-Fermion interaction (195) in the Hartree approxi- 
mation. Then the only velocity-dependent contributions are due to the current 
couplings mediated by uj and p exchanges. The corresponding contributions 
to the Landau parameters are 



and 



m = -Cl% (198) 



The derivations of relativistic F\{lo) and F[{p) are shown in Appendix C. 

Now the calculation of the convection current and the nucleon effective mass 
with the interaction (195) in the Landau method goes through the same way 
as in the case of the pion. The net result is just Eq. (185) including the 
contribution of the contact interactions (197,198), i.e., 



Fi = Fi(7r)+Fi(w), (199) 

F[ = F{(n) + F[(p). (200) 

Similarly, the nucleon effective mass is determined by (41) with 

Fi = F!(7r) + Fi(w). (201) 



In chiral Lagrangina approach the most efficient way to bring in the vector 
mesons into the chiral Lagrangian is to implement BR scaling in the param- 
eters of the Lagrangian. We shall take the masses of the relevant degrees of 
freedom to scale in the manner of BR as (128). Note again that M* is a BR 
scaling nucleon mass which will turn out to be different from the Landau ef- 
fective mass m* L [21]. For our purpose, it is more convenient to integrate out 
the vector and scalar fields and employ the resulting four-Fermi interactions 
(195). The coupling coefficients are modified compared to Eq. (196), because 
the meson masses are replaced by effective ones: 



i 



C 2 M = with M — <f>, p, uj. (202) 
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The coupling constants may also scale [22] but we omit their density depen- 
dence for the moment. 

The first thing we need is the relation between the BR scaling factor $ which 
was proposed in [14] to reflect the quark condensate in the presence of matter 
and the contribution to the Landau parameter iq from the isoscalar vector 
(uj) meson. For this we first calculate the Landau effective mass m* N in the 
presence of the pion and u fields [21] 

^ = 1 + \{FM + F ± (n)) = (l - + iMvr)))" 1 . (203) 

Next we compute the nucleon self-energy using the chiral Lagrangian. Given 
the single quasiparticle energy 

z P = jp; + c ^ N } + ^(p), (204) 

we get the effective mass as in [21] 



m L 



rriN \ dp 



p=k F j 



(205) 



from Eqs. (128), (183), and (191). Comparing (203) and (205), we obtain the 
important result 



FiH = 3(1 - (206) 

This is an intriguing relation. It shows that the BR factor, which was originally 
proposed as a precursor manifestation of the chiral phase transition charac- 
terized by the vanishing of the quark condensate at the critical point [14], is 
intimately related (at least up to p p ) to the Landau parameter iq, which 
describes the quasiparticle interaction in a particular channel. We believe that 
the BR factor can be computed by QCD sum rule methods or obtained from 
current algebra relations such as the GMOR relation evaluated in medium. As 
was shown in [21], Eq. (206) implies that the BR factor governs in some, per- 
haps, intricate way low-energy nuclear dynamics. The equivalence discussed 
above between the physics of the vacuum property $ and that of the quasipar- 
ticle interaction F 1 due to the massive vector-meson degree of freedom suggests 
that the "bottom-up" approach - going up in density with a Lagrangian whose 
parameters are fixed at zero density - and the "top-down" approach - extrap- 
olating with a Lagrangian whose parameters are fixed at some high density 
- can be made equivalent at some intermediate point. If this is so in the hot 
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and dense regime probed by relativistic heavy ion collisions, then the CERES 
data should also be understandable in terms of hadronic interactions without 
making reference to QCD variables. Because of the complexity of hadronic 
descriptions, it will be difficult to relate the two directly but the recent al- 
ternative explanation of the CERES data in terms of "melting of the vector 
mesons" inside nuclear matter manifested in the increased width of the mesons 
due to hadronic interactions [82] may be an indication for a possible "dual" 
description at low density between what is given in QCD variables (e.g., quark 
condensates) and what is given in hadronic variables (e.g., the Landau param- 
eter), somewhat reminiscent of the quark-hadron duality in heavy-light-quark 
systems [93] . A possible mechanism that could make the link between the two 
descriptions was suggested recently by Brown et al. [65] and by Kim et al. [66]. 

In the presence of the BR scaling, a non-interacting nucleon in the chiral 
Lagrangian propagates with a mass M*, not the free-space mass M. Thus, 
the single-particle current Fig. 10a is not given by (193) but instead by 



Ji-body - M „ 2 • (207) 

Now the current (207) on its own does not carry conserved charge as long 
as M* ^ M. This means that two-body currents are indispensable to restore 
charge conservation. Note that the situation is quite different from the case 
of Fermi-liquid theory. In the latter case, the quasiparticle propagates with 
the Landau effective mass m* L and it is the gauge invariance that restores m* L 
to M. In condensed matter physics, this is related to a phenomenon that the 
cyclotron frequency depends on the bare mass, not on Landau effective mass. It 
may be referred to as Kohn effect [94]. The bare mass in Kohn effect is restored 
due to the quasiparticle interactions with Galilean invariance in the same way 
as for the convection current in Section 6.1.1 [95]. This clearly indicates that 
gauge invariance is more intricate when BR scaling is implemented. Indeed 
if the notion of BR scaling and the associated chiral Lagrangian are to make 
sense, we have to recover the charge conservation from higher-order terms in 
the chiral Lagrangian. This constitutes a strong constraint on the theory. 

Let us now calculate the contributions from the pion and heavy meson degrees 
of freedom. The pion contributes in the same way as before, so we can carry 
over the previous result of Fig. 10b, 



JUody = Jfil&W ~ A(7T))7S. (208) 

This is of the same form as (194) obtained in the absence of BR scaling. It 
is in fact identical to (194) if we assume that one-pion-exchange graph does 
not scale in medium at least up to nuclear matter density. This assumption 
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Fig. 11. (a) Feynman diagram contributing to the EM convection current from 
four-Fermi interactions corresponding to the u) and p channel (contact interaction 
indicated by the blob) in effective chiral Lagrangian field theory. Th N denotes the 
anti-nucleon state that is given in the chiral Lagrangian as a 1/M correction and the 
one without arrow is a Pauli-blocked or occupied state, (b) The equivalent graph in 
heavy- fermion formalism with the anti-nucleon line shrunk to a point. 

is supported by observations in pion-induced processes in heavy nuclei. This 
means that the observation that the one-pion-exchange potential does not 
scale implies that the constant g\j '/* remains unsealing at least up to normal 
nuclear matter density with non-scaling pion mass. In what follows, we will 
make this assumption implicitly. 

The contributions from the vector meson degrees of freedom are a bit trickier. 
They are given by Fig. 11. Both the u (isoscalar) and p (isovector) channels 
contribute through the antiparticle intermediate state as shown in Fig. 11a. 
The antiparticle is explicitly indicated in the figure. However in the heavy- 
fermion formalism, the backward-going anti-nucleon line should be shrunk 
to a point as Fig. lib, leaving behind an explicit 1/M dependence folded 
with a factor of nuclear density signaling the 1/M correction in the chiral 
expansion. One can interpret Fig. 11a as saturating the corresponding counter 
term although this has to be yet verified by writing the full set of counter terms 
at the same order. These terms with Fig. 11a 



where Fi(u) and F[(p) are given by Eqs. (197,198) with M replaced by M*. 
Their relativistic forms are given in Appendix C. The total current given by 
the sum of (207), (208), (209) and (210) precisely agrees with the Fermi-liquid 
theory result (185) when we take 




k 1 ~ 



(209) 



(210) 




(211) 
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(a) 



(b) 



Fig. 12. Particle-hole contributions to the convection current. Here backward-going 
nucleon line iV -1 denotes a hole. These graphs vanish in the q/u —* limit. 



The way in which this precise agreement comes about is nontrivial. What hap- 
pens is that part of the u channel restores the BR-scaled mass M* back to the 
free-space mass M in the isoscalar current. (It has been known since sometime 
that something similar happens in the standard Walecka model (without pions 
and BR scaling) [96]). Thus, the leading single-particle operator combines with 
the sub-leading four-Fermi interaction to restore the charge conservation as 
required by the Ward identity. This is essentially the "back-flow mechanism" 
which is an important ingredient in Fermi-liquid theory. We describe below 
the standard back-flow mechanism as given in Sec. 2.1, adapted to nuclear 
systems with isospin degrees of freedom, and elucidate the connection to the 
results obtained with the chiral Lagrangian in this subsection. 

The current so constructed is valid for a process occurring very near the Fermi 
surface corresponding to the limit (u, q) — > (0, 0) where q is the spatial mo- 
mentum transfer and u is the energy transfer. In the diagrams considered so 
far (Fig. 10 and 11) the order of the limiting processes does not matter. How- 
ever, the particle-hole contribution, which we illustrate in Fig. 12 with the pion 
contribution, does depend on the order in which q = \q\ and uo approach zero. 
Thus, in the limit q/uo — > 0, the particle-hole contributions vanish whereas in 
the opposite case u/q — > 0, they do not. This can be seen by examining the 
particle-hole propagator 



g + efc - efc+g + i5 qo + e fe - e k+q - i5 

where (go, q) is the four-momentum of the external (EM) field. This vanishes 
if we set g — > with g non-zero but its real part is non-zero if we interchange 
the limiting process since for g = we have 




(212) 





(213) 
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^pf*'-')' (214) 

Figure 12 was computed by several authors (e.g., see [88]) and is given in the 
limit u/q — > by 



<V = " EM 1 ) • ^r(2)> / ^P^f - \p\)fl (215) 

where /J is the spatial and spin part of the quasiparticle interaction which is 
(188) and (B.7) without isospin part (3 — r • r')/2. The isospin factor is given 
by the Fierz transformation: 



r' " t' 

= |-^. (216) 

Note that the factor | comes from f w and |r 3 from f^. In the limit that we 
are concerned with (i.e. T = and u/q — > 0), the particle- hole contribution 
to the current is 



J P h = -^- 2 kkUfi + f[r 3 ) 

k ( F 1 (tt) + F{(7t)t 3 
M \ 6 

This holds in general regardless of what is being exchanged as long as the 
exchanged particle has the right quantum numbers. Contributions from heavy- 
meson exchanges are calculated in a similar way. Adding the particle-hole 
contribution (217) to the Fermi-liquid result (185) we obtain the current of a 
dressed or localized quasiparticle 



(217) 



= (218) 

Note that J p h precisely cancels 5 J, Eq. (180). The current Ji oc qp is the to- 
tal current carried by the wave packet of a localized quasiparticle with group 
velocity \ F = Mr- However, the physical situation corresponds to homoge- 
neous (plane wave) quasiparticle excitations. The current carried by a lo- 
calized quasiparticle equals that of a homogeneous quasiparticle excitation 
modified by the so called back-flow current [97]. The back-flow contribution 
(JiocQP — J lm) is just the particle-hole polarization current in the u/q — > 
limit, Eq. (217). 
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6.1.4 Phenomenological test 



It is not obvious that the effective nucleoli mass computed in the chiral La- 
grangian approach is directly connected to a measurable quantity although 
quasielastic electron scattering from nuclei does probe some kind of effective 
nucleon mass and Walecka model describes such a process in terms of an ef- 
fective mass. To the extent that the bulk of is related to the condensate 
through BR scaling as we can see in (205), the effective mass in the chiral 
Lagrangian can be related to the quantity calculated in the QCD sum rule 
approach for in-medium hadron masses. In BR scaling, the parameter <£> is re- 
lated to the scaling of the vector meson (p) mass. There are several QCD sum 
rule calculations for the p meson in-medium mass starting with [98]. The most 
recent one which closely agrees with the GMOR formula in medium for the 
pion decay constant /* (see Eq. (112)) is the one by Jin and Leinweber [73]: 



We shall take this value in what follows but one should be aware of the pos- 
sibility that this value is not quite firm. A caveat to this result was recently 
discussed by Klingl, Kaiser, and Weise [71] who show that the QCD sum rule 
can be saturated without the mass shifting downward by increasing the vector 
meson width in medium. For a discussion of the empirical constraints on the 
in-medium widths of vector mesons, see Friman [99] . 

Given this, we can compute m* N using (205) for nuclear matter density since 
the pionic contribution Fi(ir) is known. One finds [21] 



This can be tested in an indirect way by looking at certain magnetic response 
functions of nuclei as described below. An additional evidence comes from 
QCD sum rule calculations. Again there are caveats in the QCD sum rule 
calculation for the nucleon mass even in free-space and certainly more so in 
medium. Nevertheless the most recent result by Furnstahl, Jin, and Leinwe- 
ber [78] is rather close to the prediction (220): 



$(p ) = 0.78 ±0.08. 



(219) 



M 



(p = Po) ~ 0.70. 



(220) 




(221) 



If one writes the gyromagnetic ratio gi as 



80 



9i = 



I + 7-3 
2 



+ S91 



(222) 



then the chiral Lagrangian prediction is 



5g l = hF[-F 1 )r 3 = ^h- 1 - 



1 " ^iWj r ->- 



(223) 



In writing the second equality we have used (189), (206) and the nonet relation 
F'{p) = F(u)/9. At nuclear matter density, we get, using (219), 



Sgi(po) w 0.23r 3 . 



(224) 



This agrees with the value extracted from the dipole sum rule in 209 Bi [100], 



and agrees roughly with magnetic moment data in heavy nuclei. Nuclear mag- 
netic moments are complicated due to conventional nuclear effects. To make a 
meaningful comparison, one would have to extract all "trivial" nuclear effects 
and this operation brings in inestimable uncertainties. It should be stressed 
that the gyromagnetic ratio provides a test for the scaling nucleon mass at 
p ~ Pq. It also gives a check of the relation between the baryon property 
on the left-hand side of Eq. (206) and the meson property on the right-hand 
side. Instead of using (219) as an input to calculate 5gi, we could take the 
experimental value (225) to determine, using (223), the BR scaling factor $ 
at p pq. We would of course get (219), a value which is consistent with 
what one obtains in the QCD sum rule calculation and also in the in-medium 
GMOR relation. 

Recall that because of the pions which provide (perturbative) non-local in- 
teractions to the Landau interaction, the Landau mass for the nucleon scales 
differently from that of the vector mesons. (See (112) and (205)). This dif- 
ference is manifested in the skyrmion description by the fact that the coeffi- 
cient of the Skyrme quartic term must also scale. In the original discussion 
of the scaling based on the quark condensates using the trace anomaly [14], 
the Skyrme quartic term was scale-invariant and hence the corresponding g\ 
was non-scaling. So the scaling implied by (205) indicates that the scaling of 
g\ is associated with the pionic degrees of freedom. This is consistent with 
the description based on the Landau-Migdal g' Q interaction between a nucleon 
and a A resonance [101-103] and also with the QCD sum rule description of 
Drukarev and Levin [104] who attribute about 80% of the quenching to the 
A - N effect. 

If we equate the skyrmion relation [21,74] 



iroton 



0.23 ±0.03 



(225) 
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m* 



N 




M 

to (205), we get 



(226) 



!=(l + ^W) =(l-!fi(,r)*)" . (227) 
At nuclear matter density, this predicts 

9a(Po) » 1 (228) 

and 

iA ~ Ll = $. (229) 

We will use this relation in deriving (254). It should be emphasized that this 
relation, being unrelated to the vacuum property, cannot hold beyond p ~ po- 
Indeed as suggested by the scaling given in [14], g\(p) ~ olqa with a a constant 
independent of density, for p > p . It would be a good approximation to set 
g\ equal to 1 for p > p . 

Since one expects that when chiral symmetry is restored, qa will approach 1, 
it may be thought that the evidence for g\ « 1 in nuclei is directly connected 
with chiral restoration. This is not really the case. Neither in the skyrmion 
picture nor in QCD sum rules is the quenching of qa simply related to a 
precursor behavior of chiral restoration. This does not however mean that the 
quenching of qa carries no information on the chiral symmetry restoration. 
As suggested recently by Chanfray, Delorme and Ericson [105], if one were to 
compute all pion-exchange-current graphs at one-loop order that contribute to 
the in-medium gA, the effect of medium-induced change in the quark condesate 
would be largely accounted for. In a way, this argument is akin to that for the 
Cheshire-Cat (or dual) phenomenon we are advocating in the description of the 
quark condensate in terms of quasiparticles. Another issue that has generated 
lots of debate in the past and yet remains confusing is the interpretation 
of an effective constant g e A ~ 1 actually observed in medium and heavy 
nuclei. The debate has been whether the observed "quenching" is due to "core 
polarizations" or "A- hole effect" (or other non-standard mechanisms). Our 
view is that in the presence of BR scaling, both are involved. In light nuclei 
in which the Gamow- Teller transition takes place in low density, the tensor 
force is mainly operative and the core polarization (i.e., multiparticle-multihole 
configurations) mediated by this tensor force is expected to do most of the 
quenching, while the A-hole effect directly proportional to density is largely 
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(a) 



(b) 



Fig. 13. Examples of the second order core polarization contribution to 
Gamow- Teller transition. Downward nucleon line denote a hole. (a)2plh (b)3p2h 




Fig. 14. The resonance-exchange graph in [103] for four-Fermi contact interaction 
contribution. Its Hartree contribution decreases axial vector coupling constant in 
medium. 



P, S 7T 



A, N* 



(a) 



(&) 



Fig. 15. The resonance-exchange graphs in [103] for one-pion-exchange contribution. 
Their Fock contributions enhance the axial vector coupling constant with incorpo- 
rating the short-range correlation between nucleons. 

suppressed. The typical diagrammatic representations for the second order 
core polarization is shown in Fig. 13. In heavy nuclei, on the other hand, the 
tensor force is quenched due to BR scaling, rendering the core polarization 
mechanism ineffective while the increased density makes the A-hole effect 
dominant. Recently Park, Jung, and Min [103] applied chiral perturbation 
theory to calculate g\ at normal nuclear matter density po- The resonance- 
exchange graphs that contribute are shown in Figs. 14 and 15 and the Landau- 
Migdal g' Q effect contains both. The Hartree contribution from Fig. 14, i.e. 
A-hole effect, makes g\ quenched and the Fock contribution from Fig. 15 
enhances g\. The magnitude of quenching is two or three times larger than 
the enhancement. What is seen in nature, in our view, is the interplay between 
these two. 



The second form of (227) shows that the quenching of gA in matter is quite 
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complex, both the pionic effect and the vacuum condensate effect being con- 
founded together. Again for the reason given above, this relation cannot be 
extended beyond the regime with p p . We have no understanding of how 
this formula and the A-hole mechanism of [101,102] are related. Our effort thus 
far has met with no success. Understanding the connection would presumably 
require the short- distance physics implied by both the Landau-Migdal g' Q in- 
teraction and the Skyrme quartic term (which is known to be more than just 
what results when the p meson is integrated out of the chiral Lagrangian). 

6.2 Axial charge transition 

No one has yet derived the analogue to (185) for the axial current. Attempts 
using axial Ward identities in analogy to the EM case have not met with 
success [106]. The difficulty has presumably to do with the role of the Gold- 
stone bosons in nuclear matter which is not well understood. In this subsec- 
tion, we analyze the expression for the axial charge operator obtained by a 
straightforward application of the Fermi-liquid theory arguments of Landau 
and Migdal and compare this expression with that obtained directly from 
the chiral Lagrangian using current algebra. For the vector current we found 
precise agreement between the two approaches. 

6.2.1 Applying Landau quasiparticle argument 

The obvious thing to do is to simply mimic the steps used for the vector 
current to deduce a Landau-Migdal expression for the axial charge operator. 
We use both methods developed above and find that they give the same result. 

In free space, the axial charge operator nonrelativistically is ~ er ■ v where 
v = k/M is the velocity. In the infinite momentum frame, it is the relativistic 
invariant helicity a ■ v. It is thus tempting to assume that near the Fermi 
surface, the axial charge operator for a local quasiparticle in a wave packet 
moving with the group velocity yp = k/m* L is simply ~ a ■ vp. This suggests 
that we take the axial charge operator for a localized quasiparticle to have the 
form 

A ^o cQ p = 9a^— T -. (230) 
m L z 

As in the vector current case, we take (230) to be the oo/q — > limit of the axial 
charge operator. The next step is to compute the particle-hole contribution to 
Fig. 12 (with the vector current replaced by the axial current) in the u/q — > 
limit. A simple calculation gives 
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K. = (231) 

with 



A ' = ^0^" - h) (232) 

where I\ was defined in (190) and 



I = [ dx = In ( 1 + ^) . (233) 

In an exact parallel to the procedure used for the vector current, we take the 
difference 



Aoi ocQ P - A ; h (234) 
and identify it with the corresponding "Landau axial charge" (LAC): 



A olac = 9a^^(1 + A>). (235) 



Let us now rederive (235) with an argument analogous to that proven to be 
powerful for the convection current. We shall do the calculation using the pion 
exchange only but the argument goes through when the contact interaction 
(195) is included. We begin by assuming that the axial charge - in analogy to 
(172) for the convection current - takes the form, 



Ao z =9aJ2 J CT ' ( V P e P )n p j (236) 
where n p and e p are 2x2 matrices with matrix elements 



[n p (r, tj\ aa , = n p (r, t)5 aa < + s p (r, t) • cr aa ,, (237) 

and 

[e„(r, 0W = £ P (r, t)5 aa > + r) p (r, t) ■ v aa , (238) 

with 
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S p( r > *) = ^ 2 °"aa' *)W ( 239 ) 

In general n = 4 in the spin-isospin space. But without loss of generality, we 
could confine ourselves to n = 2 in the spin space with the isospin operator 
explicited as in Eq.(236). Then upon linearizing, we obtain from (236) 



V = <?a£ / (7^( CT - <y£) Sn *r - a- ■ (V p n° p )5e paT ) T - + ■ ■ ■ (240) 



where 



1,3, ^l + ^sr* 



H, CTT = -5 3 (p-fc)^-^- (241) 



and 



/dp 1 
To \^ fpvT,p'c>T>5np'c'T'- (242) 

■ . V ^ ) 

in analogy with (178). The equation (240) is justified if the density of polarized 
spins is much less than the total density of particles (assumed to hold here). 
The first term of (240) with (241) yields the quasiparticle charge operator 



while the second term represents the polarization of the medium, due to the 
pion-exchange interaction (188) 



5A ^ = 9A ^ T -A'. (244) 

The sum of (243) and (244) agrees precisely with the Landau charge (235). 

It is not difficult to take into account the full Landau-Migdal interactions (2) 
which includes the one-pion-exchange interaction as well as other contributions 
to the quasiparticle interaction. Thus, the general expression is obtained by 
making the replacement 



1 10 4 2 

A' -> -G' - —H' + -Hi - —H' (245) 
3 1 3 3 1 15 2 v ; 
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in (244). This combination of Fermi-liquid parameters corresponds to a I = 
£' = 1 , J = distortion of the Fermi sea [34] . We will see later that the result 
obtained with the naive Landau argument may not be the whole story, since 
the one-pion-exchange contribution disagrees, though by a small amount, with 
the chiral Lagrangian prediction derived below. 



6.2.2 Chiral Lagrangian prediction 

We now calculate the axial charge using our chiral Lagrangian that reproduced 
the Landau-Migdal formula for the convection current. Consider first only the 
pion-exchange contribution. In this case we can take the unperturbed nucleon 
propagator to carry the free space mass M. The single-particle transition 
operator corresponding to Fig. 10a is given by 



Ao\- bod y=9A^j r ^. (246) 

There is no contribution of the type of Fig. 10b because of the (G-)parity 
conservation. The only contribution to the two-body current comes from Fig. 
10c and is of the form [107] 



A-02-body — 9 A 



<7 ■ k T % 

M ~2 



with 



(247) 



2 

The factor (l/g A ) in (248) arose from replacing -4- by using the Goldberger- 
Treiman relation. 

Now consider what happens when the vector degrees of freedom are taken into 
account. Within the approximation adopted, the only thing that needs be done 
is to implement the BR scaling. The direct intervention of the vector mesons p 
and uj in the axial-charge operator is suppressed by the chiral counting, so they 
will be ignored here. This means that in the single-particle charge operator, 
all that one has to do is to replace M by M* = M$ in (246): 



body 



9 A 



a ■ kr l 



(249) 
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and that in the two-body charge operator (247), f n should be replaced by f n <& 
and M by M$: 



Aol- body = 9a^ T -A. (250) 

In the two-body operator, there is a factor (qa/ f-n) coming from the ttNN 
vertex which as mentioned before, is assumed to be non-scaling at least up 
to nuclear matter density [108], in consistency with the observation that the 
pion-exchange current does not scale in medium. 

The total predicted by the chiral Lagrangian (modulo higher-order corrections) 
is then 



/ i 

g A ^—— — (l + A). (251) 
which differs from the charge operator obtained by the Landau method, (235). 



6.2.3 Comparison between vector and axial current 

An immediate question (to which we have no convincing answer) is whether 
or not the difference between the two approaches - the Fermi-liquid vs. the 
chiral Lagrangian - is genuine or a defect in either or both of the approaches. 
One possible cause of the difference could be that both the assumed localized 
quasiparticle charge, Eq.(230), and the effective axial charge, Eq.(236), are 
incomplete. We have looked for possible additional terms that could contribute 
but we have been unable to find them. So while not ruling out this possibility, 
we turn to the possibility that the difference is genuine. 

It is a well-known fact that the conservation of the vector current assures that 
the EM charge or the weak vector charge is gy = 1 but the conservation of the 
axial vector charge does not constrain the value of the axial charge qa , that 
is, g A can be anything. This is because the axial symmetry is spontaneously 
broken. In the Wigner phase in which the axial symmetry would be restored, 
one would expect that g\ — 1- It therefore seems that the Goldstone structure 
of the "vacuum" of the nuclear matter is at the origin of the difference. 

To see whether there can be basic differences, let us look at the effect of 
the pion field. The cancellation between the two-body current J^-body (208) 
and Jp h (217) leaving only a term that changes M* to m* L in the one-body 
operator with a BR scaling mass, Eq.(207), in the EM case can be understood 
as follows. Both terms involve the two-body interaction mediated by a pion- 
exchange. It is obvious how this is so in the latter. To see it in the former, 
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we note that it involves the insertion of an EM current in the propagator of 
the pion. Thus the sum of the two terms corresponds to the insertion of an 
EM current in all internal hadronic lines of the one-pion exchange self-energy 
graph of the nucleon. The two-body pionic current - that together with the 
single-particle current preserves gauge invariance - is in turn related to the 
one-pion-exchange potential V n . Therefore what is calculated is essentially an 
effect of a nuclear force. Now the point is that the density-dependent part 
of the sum (that is, the ones containing one hole line) - apart from a term 
that changes M* to m\ in (207) - vanishes in the u/q — » limit. In contrast, 
the cancellation between (244) and (231) in the case of the axial charge, has 
no corresponding interpretation. While the one-pion-exchange interaction is 
involved in the particle-hole term (231), (244) cannot be interpreted as an 
insertion of the axial vector current into the pion propagator since such an 
insertion is forbidden by parity. In other words, Eq.(244) does not have a 
corresponding Feynman graph which can be linked to a potential. We interpret 
this as indicating that there is no corresponding Landau formula for the axial 
charge in the same sense as in the vector current case. 

In a chiral Lagrangian formalism, each term is associated with a Feynman 
diagram. As mentioned, there is no contribution to the convection current from 
a diagram of the type Fig. 10c (apart from a gauge non-invariant off-shell term 
which cancels the counter part in Fig. 10b). Instead this diagram renormalizes 
the spin gyromagnetic ratio. In contrast, the corresponding diagram for the 
axial current does contribute to the axial charge (247). As first shown in 
[89], the contribution from Fig. 10c for both the vector current and the axial- 
vector current is current algebra in origin and constrained by chiral symmetry. 
Furthermore it does not have a simple connection to nuclear force. While the 
convection current is completely constrained by gauge invariance of the EM 
field, and hence chiral invariance has little to say, both the EM spin current 
and the axial charge are principally dictated by the chiral symmetry. This 
again suggests that the Landau approach to the axial charge cannot give 
the complete answer even at the level of quasiparticle description. There is 
however a caveat here: in the Landau approach, the nonlocal pionic and local 
four-Fermion interactions (195) enter together in an intricate way as we saw 
in the EM case. Perhaps this is also the case in the axial charge, with an 
added subtlety due to the presence of Goldstone pions. It is possible that the 
difference is due to the contribution of the four-Fermion interaction term to 
(245) which cancels out in the limit u/q — > but contributes in the q/uj — > 
limit. This term cannot be given a simple interpretation in terms of chiral 
Lagrangians. Amusingly the difference between the results (see below) turns 
out to be small. 



89 



6.2.4 Numerical comparison 



To compare the two results, we rewrite the sum of (243) and (244), i.e., "Lan- 
dau axial charge" (LAC), using (41) and (189) 



A)W = ^^|f(l + A) (252) 



where 



and the sum of (249) and (250), i.e., the "current-algebra axial charge" (CAAC), 
as 



AohAAc = 9A^j(l + A) (254) 



where 



We shall compare A and A for two densities p = |po (&f = 1.50m„-) and p = po 
(hp = 1.89m,,-) where po is the normal nuclear matter density 0.16/fm 3 . 

For numerical estimates, we take 



$(p) = (l + 0.28— ^ (256) 
V Po) 

which gives $(po) = 0.78 found in QCD sum rule calculations [22]. Somewhat 
surprisingly, the resulting values for A and A are close; they agree within 10%. 
For instance at p ps Po/2, A ps 0.48 while A ps 0.43 and at p ~ po> A ps 0.56 
while A ^ 0.61. Whether this close agreement is coincidental or has a deep 
origin is not known. 



6.2.5 Test: axial charge transition in heavy nuclei 
The axial charge transition in heavy nuclei 

A{J + ) «-> B(J-) (257) 
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with change of one unit of isospin AT = 1 provides a test of the axial charge 
operator (254) or (252). To check this, consider the Warburton ratio e M EC [109] 



€mec = M exp /M sp (258) 

where M exp is the measured matrix element for the axial charge transition 
and M sp is the theoretical single-particle matrix element. There are theoretical 
uncertainties in defining the latter, so the ratio is not an unambiguous object 
but what is significant is Warburton's observation that in heavy nuclei, tMEC 
can be as large as 2: 



e HeavyNuclei = 1Q _ 2.0. (259) 

More recent measurements - and their analyses - in different nuclei [110] 
quantitatively confirm this result of Warburton. 

To compare our theoretical prediction with the Warburton analysis, we calcu- 
late the same ratio using (254) 



e< AiEC = + A). (260) 

The formula(254) differs from what was obtained in [111] in that here the non- 
scaling in medium of the pion mass and the ratio QaI f-K is taken into account. 
We believe that the scaling used in [111] (which amounted to having A/<3> in 
place of A in (260)) is not correct. 

The enhancement corresponding to the "Landau formula" (252) is obtained 
by replacing A by A in (260). Using the value for <£> and A at nuclear matter 
density, we find 



cSboW 2.1 (2.0) (261) 

in good agreement with the experimental results of [109] and [110]. Here the 
value in parenthesis is obtained with the Landau formula (252). The difference 
between the two formulas (i.e., current algebra vs. Landau) is indeed small. 
This is a check of the scaling of f n in combination with the scaling of the 
Gamow- Teller constant qa in medium. 



7 Summary 



An attempt is made and some success is obtained in this review to relate an 
effective chiral Lagrangian to an effective field theory for nuclear matter. The 
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aim is to bridge between what we know at normal nuclear density and what 
can be expected under the extreme condition, relevant in neutron stars and 
in relativistic heavy ion collisions. Furnstahl, Serot and Tang's effective chiral 
model Lagrangian FTS1 [11], which describes successfully the phenomenology 
of finite nuclei and infinite nuclear matter, is taken to imply that an effective 
chiral Lagrangian calculated in high chiral orders corresponds to Lynn's chi- 
ral soliton with the chiral liquid structure [39] in mean field. This provides 
the ground state around which quantum fluctuations can be calculated. Note 
that FTS1 is simply one of the available theories that are consistent with the 
symmetries of QCD and successful phenomenologically. We do not imply that 
FTS1 is the best one can construct as an effective theory of nuclear matter. 

The scalar sector in FTS1 develops a large anomalous dimension, which is 
interpreted as a signal of a strong coupling situation. It is suggested that the 
strong coupling theory can be transformed into a weak coupling theory if the 
chiral Lagrangian is rewritten in terms of the parameters given by BR scaling. 
A simple model, whose mass parameters are BR-scaled, is constructed and is 
shown to describe ground state properties of nuclear matter very well with fits 
comparable to the full FTS1 theory. The simple BR-scaled Lagrangian gives 
the background at any arbitrary density around which fluctuations can be 
calculated. Tree diagrams yield the dominant contributions. It is shown that 
we can make simple Walecka-type models, including our simple model. The 
models are thermodynamically consistent and the dependence of parameters 
on density are represented by the interactions of hadrons. One can also map 
the density-dependent model Lagrangian into relativistic Landau Fermi-liquid 
theory. Thus a quasiparticle picture of a strongly correlated system at densities 
away from the normal nuclear matter density is obtained. 

The BR scaling parameter $ has been identified with a Landau Fermi-liquid 
parameter by means of nuclear responses to the EM convection current. The 
Landau effective mass of the nucleon m\ is given in terms of $ and pion 
cloud, i.e. the Goldstone boson of the broken chiral symmetry, through the 
Landau parameter F*. The relation between the exchange current correction 
to the orbital gyromagnetic ratio Sgi and m* L provides the crucial link between 
$ and F± which comes from the massive degree of freedom in the isoscalar 
vector channel dominated by the uj meson. The axial charge transition in 
heavy nuclei provides a relation between $ and the in-medium pion decay 
constant f*/ f n - These relations are found to be satisfied very accurately and 
to connect physics of relativistic heavy ion collision data, e. g., dilepton data of 
CERES and nucleon and kaon flow data of FOPI(4 pi multiparticle detector) 
and KaoS(Kaon Spectrometer), etc. to low energy spectroscopic properties, 
e.g., m* L , Sgi, etc. in heavy nuclei via BR scaling. 
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8 Open issues 

While as an exploration our results are satisfying, there are several crucial 
links that remain conjectural in the work and require a lot more work. We 
mention some issues for future studies. 

• We have not yet established in a convincing way that a nontopological 
soliton coming from a high-order effective chiral Lagrangian accurately de- 
scribes nuclear matter that we know of. The first obstacle here is that a 
realistic effective Lagrangian that contains sufficiently high-order loop cor- 
rections including non-analytic terms has not yet been constructed. Lynn's 
argument for the existence of such a soliton solution and identification with 
a drop of nuclear matter is based on a highly truncated Lagrangian (ignoring 
non-analytic terms). We are simply assuming that the FTS1 Lagrangian is 
a sufficiently realistic version (in terms of explicit vector and scalar degrees 
of freedom that are integrated out by Lynn) of Lynn's effective Lagrangian. 
To prove that this assumption is valid is an open problem. 

• We do not understand clearly the role and the origin of the anomalous 
dimension d an ~ 5/3 for the quarkonium scalar field in FTS1. It is an in- 
teresting problem how the scalar in FTS1 comes to include higher order 
interactions in its anomalous scaling dimension through decimations. And 
our argument for interpreting the FTS1 with such large anomalous dimen- 
sion as a strong-coupling theory which can be reinterpreted in terms of a 
weak-coupling theory expressed with BR scaling is heuristic at best and 
needs to be sharpened, although our results strongly indicate that it is cor- 
rect. 

• There is also the practical question as to how far in density the predictive 
power of the BR-scaled effective Lagrangian can be pushed. In our simple 
numerical calculation, we used a parameterization for the scaling function 
$(p) of the simple geometric form which can be valid, if at all, up to the 
normal matter density as seems to be supported by QCD sum rule and dy- 
namical model calculations. At higher densities, the form used has no reason 
to be accurate. By using the empirical information coming from nucleon and 
kaon flows, one could infer its structure up to, say, p ~ 3p and if our argu- 
ment for kaon condensation is correct - and hence kaon condensation takes 
place at p < 3po) then this will be good enough to make a prediction for the 
critical density for kaon condensation. In calculating compact-star proper- 
ties in supernovae explosions, however, the equation of state for densities 
considerably higher than the normal matter density, say, p > 5p is required. 
It is unlikely that this high density can be accessed within the presently em- 
ployed approximations. Not only will the structure of the scaling function $ 
be more complicated but also the correlation terms that are small perturba- 
tions at normal density may no longer be so at higher densities, as pointed 
out by Pandharipande, Pethick, and Thorsson [112] who approach the effect 
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of correlations from the high-density limit. In particular, the notion of the 
scaling function $ will have to be modified in such a way that it will become 
a non-linear function of the fields that figure in the process. This would alter 
the structure of the Lagrangian field theory. Furthermore there may be a 
phase transition (such as spontaneously broken Lorentz symmetry, Georgi 
vector limit, chiral phase transition or meson condensation) lurking nearby 
in which case the present theory would have already broken down. These 
caveats will have to be carefully examined before one can extrapolate the 
notion of BR scaling to a high-density regime as required for a reliable cal- 
culation of the compact-star structure. How the scaling parameters extrap- 
olate beyond normal nuclear matter density is not predicted by theory and 
should be deduced from lattice measurements and heavy-ion experiments 
that are to come. Corrections to BR scaling as massive mesons approach 
on-shell need be taken into account. The fit to the available CERES data 
indicates however that the extrapolation to higher density - perhaps up 
to the chiral phase transition - is at least approximately correct under the 
conditions that prevail in nucleus-nucleus collisions at SPS energies. How 
this could come about was discussed in [65]. 

• In addition the behavior of g* at p > p also deviates from our simple 
form. It is expected to drop more rapidly [26]. Indeed a recent calculation 
[113] of kaon attraction to 0(Q 2 ) in chiral perturbation theory that is highly 
constrained by the ensemble of on-shell kaon-nucleon data and that includes 
both Pauli and short-range correlations for many-body effects is found to 
give at most about 120 MeV attraction at nuclear matter density. Thus the 
crucial input here is the strength of the fT~-nuclear interaction in dense 
medium. The attraction decreases from the analysis of the i^-mesic atom 
by Friedman, Gal, and Batty [114] indicating the 200 MeV attraction. If the 
attraction came down to 100 ~ 120 MeV as found in [113], this would give a 
strong constraint on the constants that enter in the four-Fermi interactions 
in the chiral Lagrangian. This would presumably account for the need for a 
dropping vector coupling g* required for p > p . Moreover Kim and Lee [115] 
found recently by renormalization group analysis that the coupling constant 
gpNN drops as density increases. This crucial information is also expected 
to come from on-going heavy-ion experiments. 

• Although we show that BR scaling parameters can be written in terms of 
Landau parameter via nonrelativistic EM current, it remains to formulate 
the relativistic mapping along the line developed in Section 5.5 where ther- 
modynamic properties of a simple BR-scaled chiral model Lagrangian in 
the mean field were shown to be consistent with relativistic Landau formula 
derived in Section 2.3. This work is needed to go to higher density region. 
Such a work is in progress. Furthermore quasiquarks must become relevant 
degrees of freedom at high density. Thus quasihadron liquid is shifted to 
quasiquark liquid as density increases. The investigation of the change in 
the shift may give a good way to connect the low density physics to the 
higher one. 
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• As seen in Section 6.2 it is not clear how we deal with Goldstone bosons 
in the scheme of Landau-Migdal approach. The Landau axial charge and 
current algebra axial charge are not the same, though they give similar 
numerical values. We do not know even whether it is possible or not that 
the Landau-Migdal approach can treat Goldstone boson properly. The study 
of it will show the way to treat the Landau Fermi-liquid theory and its scope. 

• Finally one could ask more theoretical questions as to in what way our 
effective Lagrangian approach is connected to standard chiral effective the- 
ory, which does not concern scale symmetry and its anomaly, proper and if 
the theory is to be fully predictive, how one can proceed to calculate the 
corrections to the tree-diagram results we have obtained. The first issue, a 
rigorous derivation of BR scaling starting from an effective chiral action via 
multiple scale decimations required for the problem is yet to be formulated 
but the main ingredients, both theoretical and phenomenological, seem to 
be available. The second issue is of course closely tied with what the appro- 
priate expansion parameter is in the theory. These matters are addressed in 
the paper but they are somewhat scattered all over the place and it might 
be helpful to summarize them here. The answers to these questions are not 
straightforward since there are two stages of "decimation" in the construc- 
tion of our effective Lagrangian: The first is the elimination of high-energy 
degrees of freedom for the effective Lagrangian that gives rise to a soliton 
(i.e., chiral liquid) and here the relevant scale is the chiral symmetry break- 
ing scale ~ 1 GeV and the second is that given a chiral liquid which we 
argued can be identified as the Fermi-liquid fixed point, the decimation in- 
volved here is for the excitations of scale A above (and below) the Fermi 
surface for which the expansion is made in 1/N. As discussed in Section 
2.2, 1/N ~ A/kp where A is the cutoff in the Fermi system. In bringing in 
a BR-scaled chiral Lagrangian, we are relying on chiral symmetry consider- 
ations applied to a system with a density defined by nuclear matter. Thus 
the link to QCD proper of the effective theory we use for describing fluctu- 
ations around the nuclear matter ground state must be tenuous at best. As 
recently re-emphasized by Weinberg [38], low-energy effective theories need 
not be in one-to-one correspondence with a fundamental theory meaning 
that one low-energy effective theory could arise through decimation from 
several different "fundamental" theories. This applies not only to theories 
with global symmetry but also to those with local gauge symmetry. In the 
present case, this aspect is more relevant since there is a change in degrees 
of freedom between the nonperturbative regime in which we are working 
and the perturbative regime in which QCD proper is operative. 
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Appendix 

A: Effect of many-body correlations on EOS 

In this appendix, we briefly discuss the sensitivity of the EOS to the corre- 
lation parameters of (155) at a density p > po- This is shown in Fig. A.l. 
While the parameter sets Bl, B2, B3 and B4 give more or less the same equi- 
librium density and binding energy (see Table 4), the parameter set B2 has an 
instability and B4 a local minimum at ~ 2 times the normal matter density 
whereas the sets Bl and B3 give a stable state at all density, possibly up to 
meson condensations and/or chiral phase transition. It is not clear what this 
means for describing fluctuations at a density above po but it indicates that 
given data at ordinary nuclear matter density, it will not be feasible to ex- 
trapolate in a unique way to higher densities unless one has constraints from 
experimental data at the corresponding density. In our discussion, we relied 
on the data from KaoS and FOPI collaborations to avoid the fine-tuning of 
the parameters. 



B: Relativistic calculation of F[ 

In the text, the Landau parameter F[ (or f w ) was calculated nonrelativistically 
via the Fock term of Fig. B.l. Here we calculate it relativistically by Fierz- 
transforming the one-pion-exchange graph and taking the Hartree term. This 
procedure is important for implementing relativity in the connection between 
Fermi-liquid theory and chiral Lagrangian theory along the line discussed by 
Baym and Chin [36]. 

The one-pion-exchange potential in Fig. B.l is 
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Fig. A.l. E/A—M vs. p for the Bl, B2, B3 and B4 models given in Table 4 compared 
with FTS1 theory. 



K = -91nn( t 21 • T 43 ) 
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The Dirac spinors are normalized by 



(B.l) 



u\p, s)u(p, s') = 5 SS >. 
By a Fierz transformation, we have 



(B.2) 



T21 • T43 = 2(2^41^23 — T41 • T32) 



(B.3) 
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Fig. B.l. The-one-pion-exchange diagram that gives rise to F[. 



and 



M 2 7 5 MlM47 5 M3 = ^ - M 4 7 M MiM 2 7^M 3 (B.4) 

+u 4 cr /iI 'uiiZ2<VU3 + M47 m 7 5 MiM 2 7m7 5 % + W47 5 «i^27 5 «3]- 

Remembering a minus sign for the fermion exchange, we obtain the corre- 
sponding pionic contribution to the quasiparticle interaction at the Fermi 
surface, f* = -V n (p 1 = p 4 = p,p 2 = p 3 = p',p 2 = p' 2 = k 2 F ) (see (2)). 
Decomposing p as 



r = 3 T 2 T ' (fs + fv + f T + f A + f P ) (B.5) 

where S, V, T, A and P represent scalar, vector, tensor, axial vector and 
pseudoscalar channel respectively, we find 



s : 



M 4 / 2 1 
Epin 2 q 2 + m 2 



M*f 2 1 / q 2 



E 2 F m\ q 2 + ml\ 2M 2 ) 
f _ M*f 2 1 / q 2 2a' ■ pa ■ p' - a ■ pa' p-ap'a' ■ p' 

lT E 2 mlq 2 + ml\ °" 2M 2> 2M 2 

M 4 f 2 1 / , 2a pa' p' -a pa' p-a p'a' p'\ 

^ ~ E 2 F ml q 2 + ml ' ° 2M~ 2 J 

f P = 0. (B. 



with E F = \Jk 2 F + M 2 and q — \p — p'\. Thus we obtain 
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Fig. C.l. Quasiparticle interactions in vector meson channel represented by 
four-Fermi interaction. The large blob corresponds to the blob in Fig. 11. 





f_M^_ 1_ (_ __, _ q\l-CT-CT')\ 3-t-t' 

Ff g 2 + m 2 




1 /^M 2 
3 m 2 Ff, g 

In the nonrelativistic limit, Ep ~ M and we recover (188). The factor M/Ep 
comes since there is one particle in the unit volume which decreases rela- 
tivistically as the speed increases. Note that only f s and fy in (B.6) are 
spin-independent and contribute to Ff. The fs is completely canceled by the 
leading term of fy with the remainder giving Ff. In this way of deriving the 
Landau parameter F±, it is the vector channel that plays the essential role. 



C: Relativistic calculation of Fi(uo) and J^-body 



Here we compute the contribution of vector meson channel to Landau param- 
eter Fi and EM current relativistically. The way to compute the contribution 
of p- meson channel is almost the same as uj channel. So we treat here tu-meson 
channel only. 

The blob in Fig. 11a corresponding to four- Fermi interactions can be expanded 
as Fig. C.l in random phase approximation. One bubble is represented by 



IF- = lim lim2 / ^tr [fS^-fSip + q)] (C.l) 

uJ q -+0q~>0 J (27TJ 4 

where S(p) is the Fermion propagator. The factor 2 come from isospin contri- 
bution. In the presence of Fermi sea, we can divide S(p) into the four parts; 



S(p) = ^-[(^70 - P ■ 7 + M*)( 1 + ^— ) 

2uj p Pq — UJ p + to Pq — UJ p — 10 
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+ ^ + P-~f- M * } (c . 2) 
p + uj p - id 

with Up = y/M* 2 + p 2 and n p = 9(k F — p) at T = 0. The first term is the free 
particle propagator in vacuum. The second is the particle propagator for Pauli- 
blocked state in medium. The third is for hole and the fourth is for antiparticle. 
Since vacuum contribution, i.e. antiparticle-particle in vacuum contribution, 
is canceled by counter terms and particle-hole contribution vanishes in our 
limit qjojq — > 0, the antiparticle-particle in Pauli-blocked state contribution 
remains. Then (C.l) becomes 



WV = I Wf^ tT [r{0Opl ° - p ^ + M »p7o + P • 1 ~ M*) 

+ 7> p7o + p ■ 7 - M*) 7 > p7o ~ P ■ 1 + M*)] . (C.3) 

Because of rotational invariance and zero energy-momentum transfer, only U 11 
does not vanish. 



u _ 4 ) F d 3 p 3M* 2 + 2p 2 p 

n -sJ ^ ~e; (C - 4) 



with E F = yM* 2 + k 2 F . 

The quasiparticle interaction in u;- meson channel in Fig. C.l gives 



plu r~<1 

J pp' u. 



u(p)^u(p)u(p')^^u(p') - u(p)^iu(p) I £( -C 2j4*J u(p')Yu(p') 



= Cl (u(p)Yu(p)u(p')^u(p') + u(ph z u(p) i ^l^- u(p')Yu(p')^ 

= Cl-Cf-f (C.5) 
phj F 



with chemical potential p = E F + C 2 p. Thus 



And the EM current in Fig. 11a is 

cIp/e f 
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